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The  Wiener-Hopf  technique  is^^Bfirmly  established  as  a powerful  tool  for 
research  in  certain  types  of  boundary  value  problem  arising  in  acoustics. 
Typical  problems  which  may  be  solved  exactly  or  asymptotically  with  this 
technique  concern  the  sound  and  vibration  levels  generated  by  finite  or  semi- 
infinite  planar  or  cylindrical  surfaces,  of  local  or  extended  reaction, 
immersed  in  a compressible  fluid  and  subject  to  acoustic  or  mechanical  forcing. 
However,  even  the  simplest  of  these  problems  Involves  complications  which  are  .. 
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20.  ABSTRACT  (Continued) 


irrelevant  to  an  understanding  of  the  Wiener-Hopf  method  Itself  and  its 
various  extensions.  Accordingly,  this  report  was  written  in  an  attempt  to 
display  the  operation  of  the  technique  in  an  even  simpler  physical  and 
mathematical  context,  and  thereby  to  encourage  its  more  widespread  use. 

The  report  deals  with  the  application  of  Wiener-Hopf  methods  to  one- 
dimensional  wave  motions  on  strings  and  beams,  and  in  particular  with  the 
reflection  and  transmission  from  discontinuities  in  the  mechanical  proper- 
ties of  a string.  Also  Included  is  a section  illustrating  how  a generalized 
Wiener-Hopf  problem  can  be  set  up  for  a three-part  problem  involving  a 
string  of  finite  length.  Two  dimensional  wave  problems  are  then  exempli- 
fied in  a discussion  of  the  acoustic  field  generated  by  a vibrating  half- 
plane, and  the  effect  of  uniform  mean  flow  over  the  half-plane  is  Included 
to  show  how  different  types  of  "edge  condition*  may  be  accommodated^  The 
final  section  sets  out  in  detail  the  properties  of  certain  functions  arising 
very  frequently  in  application  of  Wiener-Hopf  methods  to  acoustic  problems. 
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ABSTRACT 


The  Wiener-Hopf  technique  is  now  firmly  estab- 
lished as  a powerful  tool  for  research  in  certain 
types  of  boundary  value  problem  arising  in  acous- 
tics. Typical  problems  which  may  be  solved  exact- 
ly or  asymptotically  with  this  technique  concern 
the  sound  and  vibration  levels  generated  by  finite 
or  semi-infinite  planar  or  cylindrical  surfaces, 
of  local  or  extended  reaction,  immersed  in  a com- 
pressible fluid  and  subject  to  acoustic  or  mechan- 
ical forcing.  However,  even  the  simplest  of  these 
problems  involves  complications  which  are  irrele- 
vant to  an  understanding  of  the  Wiener-Hopf  method 
itself  and  its  various  extensions.  Accordingly, 
this  report  was  written  in  an  attempt  to  display 
the  operation  of  the  technique  in  an  even  simpler 
physical  and  mathematical  context,  and  thereby  to 
encourage  its  more  widespread  use.  The  report 
deals  with  the  application  of  Wiener-Hopf  methods 
to  one-dimensional  wave  motions  on  strings  and 
beams,  and  in  particular  with  the  reflection  and 
transmission  from  discontinuities  in  the  mechan- 
ical properties  of  a string.  Also  included  is  a 
section  illustrating  how  a generalized  Wiener- 
Hopf  problem  can  be  set  up  for  a three-part  prob- 
lem involving  a string  of  finite  length.  Two- 
dimensional  wave  problems  are  then  exemplified  in 
a discussion  of  the  acoustic  field  generated  by  a 
vibrating  half-plane,  and  the  effect  of  uniform 
mean  flow  over  the  half-plane  is  included  to  show 
how  different  types  of  "edge  condition"  may  be 
accommodated.  The  final  section  sets  out  in  detail 
the  properties  of  certain  functions  arising  very 
frequently  in  application  of  Wiener-Hopf  methods 
to  acoustic  problems. 
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1.  INTRODUCTION 


The  Wlener-Hopf  technique  was  devised  in  1931  [1]  to  deal  with 
an  integral  equation  arising  in  neutron  transport  theory,  though  its 
origins — and  indeed  its  essentials — come  from  Russian  work  in  the 
1870's  on  singular  Integral  equations  [2,3].  During  the  war  the  tech- 
nique was  extensively  applied  by  Schwinger  and  his  colleagues  [4]  to 
problems  in  electromagnetic  wave  propagation,  and  much  of  the  sub- 
sequent development  of  the  method  has  taken  place  in  applications  to 
wave  diffraction  processes.  Standard  books  on  the  method  are  those 
by  Noble  [5],  Weinstein  [6],  while  several  books  (e.g..  Carrier,  Krook 
and  Pearson  [7],  Morse  and  Feshbach  [8]) have  chapters  which  attempt  to 
introduce  the  method.  In  such  introductions,  two-dimensional  boundary 
value  problems  involving  a partial  differential  equation  for  some 
field  variable  are  invariably  used  as  the  simplest  demonstration  prob- 
lems (the  Sommerfeld  problem  of  plane  wave  diffraction  by  a semi- 
infinite  rigid  screen  being  the  best  known).  Such  problems,  however, 
bring  in  at  once  a number  of  issues  which  are  Irrelevant  to  the  exposi- 
tion of  the  W-H  method;  the  introduction  of  branch  cuts  in  the  complex 
wavenumber  plane  is  one  such  issue  which — while  it  is  an  Important 
one,  and  one  which  must  be  understood  by  anyone  wishing  to  deal  with 
wave  diffraction  problems — causes  great  difficulties  for  most  students. 

Accordingly,  an  attempt  will  be  made  in  these  course  notes  to 
illustrate  the  W-H  method  in  a much  simpler  context  than  usual.  We 
shall  study  one-dimensional  time-harmonic  waves  on  strings  and  bars,  and 
in  particular  we  will  study  the  reflection  and  transmission  properties 
of  changes  in  properties  of  the  medium,  abrupt  changes  of  density  giving 
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rise  to  various  kinds  of  standard  and  generalized  forms  of  W-H  equations. 

Not  only  is  the  derivation  of  the  W-H  functional  equation  much  simpler 
than  usual  for  these  problems,  but  its  solution  is  also  much  easier,  and 
the  final  inversion  of  a Fourier  transform  Integral  can  be  readily  per- 
formed and  the  results  seen  to  correspond  with  familiar  undergraduate 
ideas. 

The  W-H  technique  is  a method  of  solving  certain  types  of  boundary 
value  problem  in  which,  typically,  we  have  information  about  the 

( 

pressure,  say, on  the  half-plane  x < 0 and  about  the  velocity  on  the  half-  j 

plane  x > 0,  and  we  cannot  solve  for  the  radiated  acoustic  field  until 
we  know  the  pressure  all  over  the  whole  boundary,  -<»  <x<  +<”,  or 
about  the  velocity  there.  In  the  W-H  method,  a single  equation  is 
derived,  relating  the  Fourier  transforms  of  the  unknown  pressure 
(for  x > 0) , of  the  unknown  velocity  (for  x < 0)  and  of  the  given 
forcing  field  arising  from  some  prescribed  force  or  source  or  incident 
field.  The  transfomeof  the  two  unknown  distributions  are  known  (partly 
because  of  information  supplied  by  the  anticipated  physical  behavior  of 
the  system  under  study)  to  have  certain  analyticity  properties  as 
functions  of  the  wavenumber  regarded  as  a complex  variable,  and  a 
certain  crucial  step  (the  W-H  method)  and  the  use  of  some  fundamental 
theorems  of  the  calculus  of  functions  of  a complex  variable  together 
enable  one  equation  to  be  solved  for  two  unknown  functions.  Then  the 
field  everywhere  can  be  found  in  terms  of  an  inverse  Fourier  integral 
which  in  many  Instances  can  be  estimated  by  stationary  phase  or  steepest 
descent/saddle  point  techniques  [e.g.,  9,10]  or  by  generalized  function 
methods  [11,12],  In  the  problems  to  be  discussed  here  such  elaborate 
methods  are  not  needed,  and  an  exact  inversion  of  the  Fourier  integrals  i 
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can  be  obtained  with  the  aid  of  elementary  residue  calculus  (see,  e.g., 
[7]). 

It  Is  usual,  and  no  doubt  more  logical,  to  start  an  introduction  to 
a technique  of  this  kind  by  first  summarizing  all  that  will  be  required 
In  the  way  of  results  from  the  theory  of  functions  of  a complex  variable. 
As  our  problems  are  extremely  simple,  both  from  the  conceptual  and  from 
the  manipulative  point  of  view,  it  seems  unnecessary  to  start  here  with 
such  a digression,  and  we  shall  introduce  the  various  Ideas  and  theorems 
as  they  arise  In  the  course  of  the  problems.  For  a rigorous  statement  of 
the  theorems  the  reader  is  referred  to  standard  references  (e.g., 
Titchmarsh  [13]).  We  should  warn  the  reader,  however,  that  the  degree 
of  rigor  which  is  usually  adopted  by  workers  in  wave  theory  is  not 
altogether  a matter  of  mathematical  pedantry.  Certain  problems.  Involving 
coupled  wave-bearing  media  In  particular,  are  extremely  delicate,  and 
require  much  more  attention  to  mathematical  rigor  than  do  the  simple 


problems  at  hand  here. 
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2.  REFLECTION  OF  WAVES  FROM  DISCONTINUITIES  ON  A STRING 

Consider  a uniform  string  of  line  density  lying  along  the 

portion  <x  < 0 of  the  x-axis  under  tension  T.  Consider  this  string 

to  be  joined  at  x = 0 by  a massless  connection  to  a string  of  density  pj 

lying  along  0 < x < +“,  the  tension  there  also  being  T.  Transverse 

1 

waves  of  arbitrary  profile  propagate  at  speed  C^  = (T/p^)^  on  the 
left-hand  string,  and  at  speed  C^  = (T/p^)^  on  the  right-hand. 
Alternatively,  with  time  dependence  exp(-liot),  lo  > 0,  understood  through- 
out, waves  in  X < 0 have  wavenumber  k = oj/C  while  those  in  x > 0 

o o 

have  wavenumber  = to/Cj.  We  wish  to  determine  the  reflected  and  trans- 
mitted waves  in  X < 0,  x > 0 respectively  when  a progressive  wave  with 

displacement  exp(ik  x)  is  incident  upon  the  junction  from  x < 0.  [See 
o 

Figure  1.] 

This  is  a trivial  problem  to  which  the  solution  can  be  found  by 
elementary  methods  and  which  can  be  generalized  to  cover  a variety  of 
different  conditions  at  the  junction.  It  and  its  generalizations  are 
also  suitable  for  introducing  the  W-H  method  very  simply.  We  start 
by  writing  the  total  displacement  as 


y + exp(ik^x) 

and  as 

y 


in  X < 0 

in  X > 0 


(2.1) 


The  reason  for  this  is  that  then  y (which  might  be  called  the  scattered 

field)  must  take  the  form  of  an  outgoing  wave  as  x +<»and  as  x -►-<», 

and  this  enables  us  to  state  something  about  the  Fourier  transform  of  y. 

We  have  first,  however,  to  make  the  wavenumbers  k and  k slightly 

o 1 


I 
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complex 


(2.2) 

where  real  and  Imaginary  parts  of  both  wavenumbers  are  both  positive 
and  where  the  Imaginary  parts  are  small  (and  in  the  end  vanishingly  small 
compared  with  the  real  parts. 

The  rationale  for  this  is  as  follows.  We  start  by  assuming  a time 
factor  exp(-i(i)t)  with  u)  > 0.  (We  could  just  as  well  take  exp(+iut),  but 
the  choice  of  exp(-l(jJt)  is  helpful  for  reasons  that  should  emerge  later.) 
Then  as  x “ the  phase  factor  of  an  outgoing  wave  must  be  exp(+ikjX) 
where  k^  = w/C^.  Giving  k^  a small  positive  imaginary  part  therefore 
makes  an  outgoing  wave  decay  as  x -►  + ®,  like  exp(-kj^x),  corresponding 
to  the  presence  of  small  internal  dissipation  in  the  string.  Similarly, 
as  X -*■  - "an  outgoing  wave  will  have  the  phase  factor  exp(-ik^x),  and 
then  will  also  be  exponentially  damped  if  we  give  k^  a small  positive 
imaginary  part. 

We  therefore  now  know  that  for  x > 0,  y(x)  is  a continuous  function 
which  decays  like  exp(-kj^x)  as  x -►  ».  Its  HALF-RANGE  FOURIER  TRANSFORM 
then  has  certain  properties  as  a function  of  complex  wavenumber  s (that 
being  a useful  symbol  for  the  wavenumber  which  does  not  have  any  particu- 
lar significance,  the  symbol  k for  example  often  being  associated  with 

a particular  wavenumber).  Define 
00 

Y^(8)  “ j y(x)  exp  isx  dx.  (2.3) 

*0 


k -*•  k + ik  , 
o or  ol 


k -►  k + ik 

“^ir  li 
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Then  Y^(s)  exists  as  an  analytic  function  of  s at  all  points  In  the 
complex  s-p’ane  for  which  the  Integral  converges.  The  Integral  over 
any  finite  range  of  x certainly  converges  (and  y(x)  Is  certainly  Inte- 
grable  near  x=0)  so  that  the  convergence  Is  dictated  by  the  behavior  of 
y(x)  as  X -►  + oo.  There 

y(x)  exp  isx  ~ exp(-(s^  + k^^)x]  (2.4) 

where  s = s^  + and  the  integral  up  to  infinity  therefore  converges 

if  s^  + > 0.  It  may  also  converge  for  some  other  values  of  s,  but 

what  can  be  guaranteed  on  the  basis  of  the  anticipated  behavior  of  y(x) 
as  X -*•  + <»  is  that 

Y^(s)  is  analytic  in  an  upper  half-plane 

R : Im  s > -k, . (2.5) 

Moreover,  Y^(s)  has  non-growing  algebraic  behavior  as  |s|  -*■  " every- 
where within  the  domain  R^  if  y(x)  is  finite  at  x = 0+  or  has  an 
integrable  singularity  at  x = CH-,  so  that  |y_^(s)1  = 0(|s|  ^)  say  for 
some  A ^ 0 as  |s|  ->•  “ along  any  radius  in  the  domain  R^.  If,  as  may 
occur  in  some  problems,  y(x)  has  a non-integrable  singularity  at  x = 0+ 

(that  is,  a singularity  at  least  as  strong  as  x~*)  then  y(x)  has  to 
be  regarded  as  a generalized  function  whose  generalized  half-range  trans- 
form Y^(s)  is  still  analytic  in  an  upper  half-plane  R^,  but  now  can  have 
algebraic  growth  as  |sl 

By  analytic  we  mean  that  Y^(s)  is  single-valued  and  has  a unique 
derivative 
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Jllm 
h -►  o 


(2.6) 


y;(s) 


Y^(s+h)  - Y^(8) 
h 


as  the  point  s+h  approaches  s along  any  path  In  the  plane. 

In  a precisely  similar  way,  the  half-range  transform  of  the  scattered 
displacement  y(x)  for  x < 0 

Y_(s)  = y(x)  exp  Isx  dx  (2.7) 

—CO 

Is  an  analytic  function  of  s In  a lower  half  of  the  s-plane, 

R ; Im  s < + k (2.8) 


and  for  a function  y(x)  Integrable  at  x = 0-  Y (s)  has  non-growing 
algebraic  behavior  at  infinity  in  R , so  that  |Y  (s) | = 0(|s|  ^)  for 
some  p >.  0 as  | s | -*■  <»  along  any  radius  in  R . (Algebraic  growth  of 
Y (s)  is  permitted  if  y(x)  has  a non-integrable  singularity  as  x 0-.) 

Since  and  k^^  are  strictly  positive,  the  FULL-RANGE  FOURIER 
TRANSFORM 


Y(s) 


Y^(s)  + Y_(s) 


00 

y(x)  exp  isx  dx 


exists  as  an  analytic  function  of  s in  a strip  D, 


(2,9) 


D - R^  n R : -k, . < Im  s < + k ^ (2.10) 

formed  by  the  intersection  of  the  overlapping  upper  and  lower  half-planes 
R^  and  R_ . [ See  Figure  2 . ] 

The  notation  and  conventions  for  the  Fourier  transforms  are  such 
that  (^functions  arise  from  half-range  trimsforms  over  the  positive 
x-axis,  ^functions  from  half-range  transforms  over  the  negative  x-axls. 
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Alternatively,  Y^(s)  can  be  regarded  as  the  full-range  transform  of  the 
generalized  function  y(x)  H(x), 


Y^(s) 


r 


00 

y(x)  H(x)  exp  isx  dx 


(2.11) 


and  correspondingly 

f+  00 

Y (s)  = I y(x)  H(-x)  exp  isx  dx  , (2.12) 

H(x)  denoting  the  Heaviside  unit  function,  H(x)  = 1 for  x > 0,  H(x)  = 0 
for  X < 0. 

The  FOURIER  INVERSION  THEOREMS  run  as  follows:-  if  is  any  contour 
from  - ® to  +°° lying  within  the  domain  R^  and  such  that  Y^(s)  exp (-isx) 
is  integrable  over  this  contour,  then 


2TT 


Y^(s) 


exp (-isx) 


ds  = y(x) 
= 0 


(x  > 0) 
(X  < 0) 


both  of  which  are  contained  in 


y(x)  H(x)  = — r Y (s)  exp(-isx)  ds  (2.13) 

2Tt  Jc^ 

The  fact  that  the  Integral  vanishes  for  x < 0 follows  from  Cauchy's 
theorem  applied  to  a closed  contour  T consisting  of  the  contour  with 
its  ends  joined  by  a large  circular  arc  in  R^.  Since  the  integrand  is 
analytic  everywhere  in  R^ 

I Y (s)  exp (-isx)  ds  * 0 

•'r 

and  the  Integral  along  the  circular  arc  vanishes  if  x < 0 
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because  exp(-l8x)  is  then  exponentially  small  at  infinity  in  and 
Y^(s)  is  at  least  algebraically  small  there.  (This  result  should  really 
be  proven  carefully  using  Jordan's  lemma;  see  [7,13].) 

Correspondingly,  if  C is  any  contour  from  -«>  to  +««>  lying  within 
R and  such  that  Y (s)  exp(-isx)  is  integrable  over  it,  then 


y(x)  H(-x) 


— f Y_(s) 

2Tr  Jc 


exp(-isx) 


ds 


(2.14) 


Now  if  C is  any  contour  from  to  +“  lying  everywhere  within  the 
strip  of  D of  overlap  between  R.  and  R , then  we  can  identify  C with  C 
for  (2.13),  C with  C for  (2.14)  and  obtain  by  addition  of  (2.13)  and  (2.14) 
the  inversion  theorem 

y(x)  * — I y(s)  exp(-lsx)  ds  (2.15) 

2Tr  Jc 


We  remark  here  that  our  convention  exp(-ia)t)  for  the  time  factor 
is  consistent  with  the  formulas  (2.9)  and  (2.13),  in  the  sense  that  we 
are  really  taking  F.T.'a  in  space  and  time  with  the  definitions 


Y(8,U)) 

y(x,t) 


y(x,t)  exp(isx  + iut)  dx  dt 


Y(8,a))  exp(-i8x  -iwt)  ds  du) 


Note  that  if  a time  factor  exp(-tlfa)t)  were  taken,  and  were  defined 

as  k^  “ then  we  would  have  to  give  k^,kj  small  negative 

Imaginary  parts  in  order  to  secure  a strip  of  overlap  in  which  to  take 
F.T. 's  in  X. 
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The  point  of  using  half-range  transforms  is  that  conditions  are 
different  according  as  x ^ 0,  so  that  a two  sided-transform  cannot 
immediately  be  applied.  For  x > 0 we  have  the  equation  of  motion 


dx^  ‘ 


1 y = 0 . 


Multiply  by  exp(isx)  and  integrate  from  0 to  Then  the  last  term 


produces 


kj  Y^(s) 


(2.16) 


provided  seR^,  while  the  first  term  can  be  integrated  by  parts  to  give 


fd^y  isx  . dy  isx  “ . . isx  ® \ 

— f e dx  = ^ e -is  [ye  -s^Y^(s) 

dx  o 

From  the  anticipated  behavior,  y~exp(  ik^^x)  as  x -► +<*>  we  see  that 

e^®*  and  ye^^^  both  vanish  as  x +®°  provided  seR  . If  we  write 
dx  + 


« » 9^  Im  « . • / /N . \ ^ inj  I ^ » 

y(o+)  = X . 04.  yM>  y = X - o+y 


(2.17) 


then  we  have 


(s^  - k|)  Y^(s)  - -y’((>+)  + isy((H) 

as  a statement,  for  seR^,  of  the  equation  of  motion  in  x > 0. 

Applying  the  same  procedure,  this  time  for  seR  , to  the  equation 


(2.18) 


k^ 

O 


y - 0 


(2.19) 


which  holds  in  x < 0, produces 


(s*  - k^)  Y (s)  - + y’(O-)  - isy(O-) 
o - 


(2.20) 
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The  system  (2.18),  (2.20)  is  completed  when  the  boundary  condition 
at  X <■  0 is  specified.  In  the  present  simplest  case,  in  which  there  can 
be  no  difference  of  transverse  force  from  x 0-  to  x ■ (H-,  the  boundary 
condition  is  that 


T — (v  + ^ Him  ^ ^ 

X -*■  0-  dx'‘^  X -*•  Of  dx  * 


(2.21) 


so  that 


y'(O-)  + ik  - y’(Of) 
o 


(2.22) 


while  the  total  displacement  of  the  string  must  also  be  continuous,  so 


Him  , . ik_x.  Him 

, - 0-  ' > ■ X - Of  y 


(2.23) 


y(O-)  + 1 - y(O)-) 


(2.24) 


Eliminate  say  y(O-)  and  y' (0-)  from  (2.20)  by  using  (2.22)  and  (2.24), 


giving 


(s^  - k*)  Y_(s)  = y'(Of)  - isy(Of)  + i(s  - k^) 


(2.25) 


and  if  we  now  add  this  is  to  (2.18),  the  unknowns  y'(0+)  and  y(0+) 
disappear  and  we  have 

(s'  - kj)  Y^(8)  + (s2  - k*^)  Y_(8)  - i(s  - k^)  . 

It  is  convenient  to  divide  through  by  (s  - k^),  which  we  may  do  because 
the  equation  Itself  is  only  meaningful  in  the  strip  D and  the  zeros  +k 

o 

lie  outside  that  strip.  Then  we  get  a standard  form  of  Wlener-Hopf 
functional  equation 
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r 


K(,)  * Y_(,)  - jjj; 


(2.26) 


where  the  kernel  Is 


K(s) 


s'  - k! 


(2.27) 


The  W-H  equation  relates  a linear  combination  of  the  unknown  half- 
range transforms  Y^(s)  to  a function  related  to  the  incident  wave  field 

exp  ik  X.  The  equation  holds  in  the  strip  of  overlap  D,  and  the  coeffi- 
o 

dents,  K(s)  and  (s  + k ) \ are  analytic  and  non-zero  in  the  interior 

o 

of  the  strip  in  the  simplest  cases,  though  cases  in  which  K(s)  has 
zeros  or  poles  in  the  interior  of  D can  also  be  handled  straighforwardly . 

In  the  next  section  we  show  how  the  general  situation  of  the  W-H 
equation  can  be  obtained  by  inspection  and  how  the  Fourier  transform 
integrals  can  be  inverted  to  yield  explicit  solutions  for  the  trans- 
mitted and  reflected  waves.  Following  that  we  look  at  differences 
which  arise  when  the  boundary  conditionsat  x = 0 are  changed,  and  when 
the  strings  are  replaced  by  elastic  beams. 


3.  SOLUTION  OF  THE  W-H  EQUATION 

Assuming  that  we  have  obtained  an  equation  (2.26)  with  the  coeffi- 
cient of  one  or  other  of  Y^(s)  reduced  to  unity,  the  first  crucial  point 
lies  in  the  W-H  FACTORIZATION  of  K(s) . In  this  we  express 


K(s)  - K_|_(s)  K_(s) 


(3.1) 


as  the  product  of  two  functions  of  which  K^(s)  is  analytic  and  non-zero 
in  and  of  at  most  algebraic  growth  at  infinity  in  R^,  while  K (s)  is 
analytic  and  non-zero  in  R and  of  at  most  algebraic  growth  at  infinity 


13 


there.  (Here  we  assume  that  K(s)  has  no  zeros  In  the  strip  D,  a point 
we  shall  examine  again  later.)  It  is  remarkable  that  such  a factoriza- 
tion exists  for  any  K(s)  which  is  analytic  and  non-zero  in  D but  which 
may  have  any  kind  of  singular  behavior  outside  D;  the  proof  was  given 
by  Wiener  and  Hopf  [1].  Here  the  truth  of  the  theorem  is  obvious; 


K^(s) 


(3.2) 


is  one  factorization  with  the  required  properties,  and  is  such  that 
K.^(s)  1 at  infinity  in  R^.  After  we  have  completed  the  solution 

we  shall  return  to  the  uniqueness  or  otherwise  of  this  factorization. 

Because  K (s)  is  free  from  zeros  in  R , and  in  particular  in  D, 
we  can  divide  (2.26)  through  by  it,  to  get 


Y (s) 

^ if® 


1 

(s  + k ) K (s) 
o - 


(3.3) 


The  analyticity  properties  of  the  terms  on  the  left  here  are  known, 

while  the  function  on  the  right  is  neither  a @ function  nor  a Qf unction. 

Our  next  object  is  to  write  it  as  the  SUM  of  the  two  functions  analytic 

in  R^  and  R respectively,  and  of  algebraic  behavior  at  " in  those 

half-planes.  Again,  the  existence  of  this  ADDITIVE  SPLIT  is  assured 

by  the  W-H  theorem  [1],  and  here  we  can  again  see  how  to  perform  the  split 

by  inspection  (and  in  this  case  the  method  of  inspection  is  widely 

useful  and  should  be  carefully  noted) . 

The  function  i/ (s  + k^)  K (s)  is  analytic  in  R_  except  for  the 

pole  at  s >•  ~k^.  Near  the  pole,  the  function  behaves  like 

1/(8  ->■  k ) K (-k  ),  so  that  we  can  isolate  the  pole  behavior  by  adding 
0 — 0 

and  subtracting  this  term,  to  give 
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(3.4) 


1 1 f 1 1 ^ , i 

(8  + k ) K (s)  * (8  + k )k  (s)  "k  (-k  )/  (s  + k^)  K (-k^) 

The  first  term  no  longer  has  a pole  at  s = -k^,  for  near  there 

^ “ (s  + k^)x  (function  analytic  at  -k^)  , 

- —o' 

and  it  is  therefore  aQfunction  G (s)  say,  while  the  correction  term 

i/(s  + k ) K (-k  ) is  evidently  a@  function,  G (s)  say.  Note  that 
0 — 0 + 

this  additive  split  is  not  restricted  to  any  particular  form  of  K.^(s), 

but  turns  on  the  presence  of  a pole  term  (s  + k )“^  only. 

o 

We  now  rewrite  the  equation  (3.3)  as 

T (s) 

K^(s)  Y^(s)  - G^(s)  - G_(s)  - , (3.5) 

and  consider  the  function  E(s)  defined  by 

E(s)  - K^(s)  Y^(s)  - C^(8)  . (3.6) 

This  is  a function  defined  and  analytic  throughout  and  of  algebraic 
behavior  (algebraic  growth  or  decay)  at  “ in  R^.  E(8)  is  not  defined 

by  (3.6)  except  in  R^.  However,  within  D,  E(s)  can  equally  well  be 
defined  by 

Y_(s) 

E(s)  - G_(a)  - (3*7) 

and  this  definition  then  CONTINUES  ANALYTICALLY  the  function  E(8), 
defined  originally  only  in  R^  by  (3.6),  through  the  strip  D of  overlap 
into  the  lower  half-plane  R , and  there  £(&)  is  also  analytic  and  of 
algebraic  growth  or  decay  at  infinity. 
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Thus  we  have  a function  E(s),  defined  by  (3.6)  In  and  (3.7)  In 
R_,  these  two  definitions  being  Identical  throughout  D,  which  Is  analytic 
and  of  at  most  algebraic  growth  at  Infinity  In  the  entire  complex  plane. 
According  to  the  extended  version  of  Liouvllle's  theorem,  the  most 
general  such  function  Is  a polynomial  P(s),  and  hence  the  most  general 
solution  to  the  W-H  equation  (2.26)  Is 


Y^.(s) 


G_^(s)  + P(s) 
K^(s) 


(3.8) 


Y (s)  - K (s)  [G  (s)  - P(s)] 


for  any  polynomial  P(s). 

In  most  applications  the  polynomial  Is  fixed  In  degree,  and  some- 
times also  completely,  from  considerations  of  "edge  conditions"  at  x *■  0. 
The  edge  conditions  determine  the  behavior  of  Y^(s)  at  infinity  in  R^. 
This  Is  seen  most  clearly  If  we  go  to  infinity  along  the  Imaginary  axis. 
Let  s > lu  with  u purely  real  and  positive.  Then 


y^(iu)  - j y(x)  exp(-ux)  dx 
o 

which  is  the  one-sided  Laplace  transform  of  y(x).  When  u ■*•  + “ only 
small  values  of  x make  any  contribution  to  the  Integral,  and  in  fact 
WATSON's  LEMMA  [see,  e.g.,9]  states  that  under  appropriate  conditions, 
the  asymptotic  behavior  of  Y^(lu)  as  u is  obtained  by  inserting 

in  (3.9)  the  asymptotic  expansion  of  y(x)  as  x 0+  and  integrating 
term  by  term.  The  process  can  also  be  used  in  an  Inverse  fashion  to 
find  the  behavior  of  y(x)  as  x (H  by  examining  the  behavior  of  Y^(lu) 
as  u +«>[14].  Note,  however,  that  the  behavior  of  y(x)  as  x -►  Oi- 
ls NOT  simply  related  to  the  behavior  of  the  full-range  transform  Y(s) 


(3.9) 
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for  large  values  of  s.  If  one  knows  Y(s),  one  has  then  to  decompose 
Y(s)  into  Y (s)  + Y (s)  and  then  look  at  Y (s)  for  large  values  of  s,  a 

+ • *T* 

procedure  explained  In  detail  In  [14]. 

In  a similar  way,  the  asymptotics  of  y(x)  as  y -*•  0-  are  determined 
by  those  of  Y (-iv)  as  v -►  + <»,  using  Watson's  lemma. 

Now  In  our  problem  we  assume  that  the  deflection  y in  the  reflected 
or  transmitted  wave  Is  finite  as  x -*■  Q±.  Then  the  leading  order  term 
in  the  expansion  of  Y^(iu)  is 

(Const)  I exp(-ux)  dx  »=  (Const)u~*  (3.10) 

o 


and  similarly  for  Y (-iv),  so  that  Y^(x)  are  each  0(s“*)  at  infinity  in 
R , respectively.  K (s)  each  tend  to  1 at  infinity,  while  G (s)  are 

+ + i 

each  0(s~M.  It  then  follows  from  (3.6)  that  E(s) , which  is  the 
polynomial  P(s),  is  0(s"M  at  infinity  in  R^,  and  from  (3.7)  that  it 
is  0(s“^)  at  infinity  in  R_,  and  hence,  because  R^  have  a common  strip, 
P(s)  is  a polynomial  which  vanishes  everywhere  at  infinity  like  s“^ . 
Therefore  P(s)  must  be  identically  zero,  and  the  solution  subject  to 
the  condition  of  flnlteness  at  the  junction  x = 0 is 


Y^(s)  - G^(s)/K_^(8) 

Y_(8)  - G_(8)  K_(s) 
and  in  explicit  form  this  gives 


(3.11) 
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21k^ 

“ (s  + k ) (k  + k ) 
1 o 1 


Y (s) 


1 

s + k 

o 


2ik 

o 

(s  + k )(k^  + k ) 

O O 1 


(3.12) 


Before  proceeding  to  the  Inversion  of  the  Fourier  transform  we 
return  to  the  question  of  the  uniqueness  of  the  K^(s).  Let  one  specific 
factorization  be  K^(s)  K (s) . Then  in  any  other  factorization,  the 
factors  can  be  written 


[A^(s)  K^(s)][A_(s)  K_(s)] 


where,  since  K^(s)  are  analytic  and  free  of  zeros  in  R^,  A^(s)  must 
also  be  analytic  and  free  of  zeros  in  R^,  respectively.  Further, 


A^(s)  A_(s)  » 1 


for  seD,  and  because  A (s)  has  no  zeros  in  R 


*+<•>  ■ TJi) 


again  for  seD.  Define 


F(s)  - A^(s)  seR^ 


F(s) 


A (s) 


seR 


» 


Then  F(s)  is  analytic  throughout  the  whole  s-plane  (an  ENTIRE  function). 
Suppose  further  that  the  factors  are  required  to  have  some  specified 
algebraic  behavior  at  infinity  In  that  respective  half-plane.  Then  A^(s) 
each  tend  to  constant  values  at  infinity  In  R^,  and  the  entire  function 
F(s)  has  constant  values  everywhere  at  infinity.  By  Liouville's  theorem 
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the  only  such  function  is  F(s)  ■ Constant  F , so  that  If  K (s) , K (s) 

o + ~ 

are  one  pair  of  factors,  any  other  pair  must  be  of  the  form 
F^  K.(s),  F-‘  K (s)  . 

In  other  words,  a factorization  with  prescribed  algebraic  behavior  at 

infinity  is  unique  up  to  multiplication  of  K^(s)  by  a constant  F^  and 

division  of  K (s)  by  the  same  constant  F . Note  the  Importance  of  the 
- o 

restriction  to  algebraic  behavior  at  infinity.  In  many  applications, 
part  of  K(s)  can  be  represented  as  an  infinite  product  from  which  the 
split  K^(s)  K_(s)  can  be  effected  by  Inspection.  Usually,  however,  the 
infinite  series  of  factors  which  gives  K^(s)  or  K (s)  has  exponential 
behavior  at  infinity  in  or  R_,  and  then  it  is  necessary  to  divide 
say  K^(s)  by  an  entire  function  with  the  appropriate  exponential 
behavior  at  infinity  in  R^,  so  that  the  resulting  factor  behaves 
algebraically,  at  the  same  time  multiplying  K (s)  by  the  same  factor 
to  eliminate  exponential  behavior  at  infinity  in  R . Several  examples 
of  this  are  given  in  the  book  by  Noble  [5], 
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4.  INVERSION  OF  THE  FOURIER  INTEGRALS 

From  (3.12),  the  full-range  transform  of  the  scattered  field  y(x) 
is  given  by 

Y(s)  * Y^(s)  + Y_(s)  , 

and 

y(x)  “ j*  exp(-isx)  ds 

C 

where  C runs  from  - «>  to  +<»in  the  strip  D.  Since  Y^(s)  behave 
algebraically  at  infinity  in  R^,  the  convergence  of  the  Integral  is 
dictated  by  the  exp(-isx)  factor. 

For  X > 0,  close  the  contour  C with  a large  semi-circle  in  the 
lower  half-plane.  The  contribution  from  this  semi-circle  vanishes 
as  the  radius  becomes  infinite  because  exp(-isx)  is  exponentially  small 
for  x > 0 and  Im  s < 0.  (Again,  a more  careful  proof  of  this  should 
really  be  given.)  Further,  inside  the  closed  contour  consisting  of  C 
and  the  large  semi-circle,  Y (s)  is  analytic  and  by  Cauchy's  theorem 
makes  zero  contribution  to  the  integral.  Y^(s)  does  have  singularities 
in  R , but  these  consist  in  fact  of  just  a simple  pole  at  s = -kj . 
Noting  that  the  contour  is  described  clockwise  rather  than  counter- 
clockwise, we  have 

y(x>0)  = (-2Tri)  (Residue  of  Y^(s)  at  s - -k^)  exp(lkjX) 

2ik 

“ (“i)  exp(ikjX) 

o 1 

“ T exp(ikjx) 

where  the  transmission  coefficient  is 
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2k 


k + k, 
o 1 


(4.1) 


For  X < 0,  complete  the  contour  C with  a large  semi-circle  in 
R^.  There  is  no  contribution  from  the  large  semi-circle,  and  applica- 
tion of  Cauchy's  theorem  now  gives 


y(x  0)  “ ^ (■*■  2TTi)  (Residue  of  Y_(s)  at  s = + k^)  exp(-  ik^x) 
(2ik  ) 

“ " (2k  Kk" -f  T)  ‘^o  - 

0 0 1 

“ R exp(-  ik  x) 
o 

with  a reflection  coefficient 
/k  - k. 


R = 


1 


k + k, 

< o 1 , 


(4.2) 


It  is  a trivial  matter  to  check  that  these  solutions  for  the  reflected 
and  transmitted  waves  satisfy  the  conditions  (2.23)  and  (2.24). 

In  the  kinds  of  problems  encountered  in  acoustics,  the  function 
Y^(s)  usually  has  a branch  point  singularity  in  R as  well  as  simple 
poles,  and  correspondingly  Y (s)  may  have  one  or  more  branch  points 
and  simple  poles  in  R^.  The  pole  contributions  are  unaffected  (except 
in  certain  critical  circumstances)  and  give  rise  to  natural  modes  of 
the  system,  analogous  to  the  reflected  and  transmitted  waves  found  here. 
When  the  dissipation  factors  are  small,  some  poles  will  lie  close  to 
the  real  axis  and  give  rise  to  propagating  modes,  as  here.  Others 
may  lie  close  to  the  imaginary  axis  and  give  non-propagating  modes,  of 
the  kind  found  in  wave-guide  problems  and  in  the  motion  of  plates  and 
beams.  In  other  problems  poles  may  be  present  in  the  complex  plane 


21 


and  yet  may  never  be  captured  In  the  appropriate  deformations  of  the 
Integration  path  which  will  have  to  be  made,  and  such  poles  then 
represent  no  distinct  and  Identifiable  structure  of  the  field.  Pole 
contributions  also  serve.  In  acoustics  problems,  to  represent  the 
abrupt  changes  that  would  occur  according  to  geometrical  optics  as  one 
crosses  boundaries  between  illuminated,  reflected,  and  shadow  wave 
zones.  These  pole  contributions  have  to  be  supplemented  In  more 
complicated  problems  by  Integrals  around  the  branch  cuts  which  join 
branch  points.  There  are  many  techniques  for  estimating  the  contributions 
from  branch  cut  integrals  [9,10,11,12]  including  cases  where  various 
kinds  of  singularities  come  close  together,  and  even  coincide.  Generally 
branch  cut  Integrals  represent  forced  near-field  behavior,  which  decays 
algebraically  away  from  a junction  or  discontinuity  of  the  two-part 
system  usually  studied  by  the  W-H  technique,  leaving  only  the  natural 
propagating  modes  at  large  distances. 
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5.  DIFFERENT  CONDITIONS  AT  x »■  0 

Suppose  now  that  the  strings  are  each  stretched  to  tension  T, 
but  are  joined  by  a mass  m which  is  free  to  slide  on  a smooth  wire 
perpendicular  to  the  strings.  The  condition  of  continuity  of  displace- 
ment of  the  string  remains  in  force,  so  that  as  in  (2.24) 


y(O-)  + 1 = y(0+)  (5.1) 

while  the  equation  of  motion  of  the  particle  is 


T ^^(0^■)  - Tt-  (y  + y(0+)  (5.2) 

dx  ax 

since  y(0+)e  is  the  particle  displacement.  We  now  find  that  it  is 
impossible  to  eliminate  all  of  y(0±),  y'(0±)from  the  equations (2. 18) , 

(2.20),  (5.1),  and  (5.2),  which  are  statements  of  the  equationsof 
motion  and  the  boundary  conditions.  The  simplest  W-H  equation  one  can 
get,  replacing  (2.26),  turns  out  to  be 


K(s)  Y^(s)  + Y_(s) 


1 mu)^ 

s + k T 
o 


y(tH-) 

s^  - 

o 


(5.3) 


containing  two  unknown  functions,  Y^(s),  and  the  unknown  constant 
y(0+)  also.  Supposing  y(0+)  were  known,  however,  we  proceed  as  before. 
The  additive  split  of  i/(s  + k ) K (s)  * G (s)  + G (s)  has  been  given 

O — T — 

in  (3.4);  for  the  other  term  in  (5.3)  we  proceed  in  a similar  fashion 
and  get 


3sli = H^(8)  + H (s)  (5.4) 

T (8  2-k2)K_(s)  + 

where 
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H^(8)  - ^ y(CH-)  (g  ^ ) (_2it  ) K (-k  ) 


(5.5) 


H (s) 


mu 


(s  + k^)  [^(s  - k^)  K_(s)  * (-2k^)  K_(-k^)j 


The  entire  function  E(s)  Is  again  zero  and  thus  we  find 

.2 


t—K  (-'k  T 


mu 


2k  T 
o 


y((H-) 


(s  + k ) K (-k  ) 
0-0 


(5.6) 


From  this  equation  we  can  now  determine  the  value  of  y(CH),  for 
we  recall  from  03  that  the  behavior  of  Y^(s)  as  s in  is  related 

to  the  behavior  of  y(x)  as  x -►  0+  — i.e.,  to  y(0+).  As  s ->■  ® in  R_j_ 
we  have 


(^(s)  - f 


y(0+)  exp(lsx)  dx  + ... 


(^)  + ^ ^ (provided  Ims  > 0) 


(5.7) 


while  K^(s)  -*■  1 as  s -►  <*>  in  R^.  Therefore  the  leading  terms  of  (5.6) 
state  that 


iy((H) 

s 


n»ii 


y(o+) 


sK  (-k  ) 2k  T s K (-k  ) 
- o o - o 


which  requires 


y(0f)  j 1 + 2k  T K (-k  ) i ” K (-k  ) 


- o 


(5.8) 


(5.9) 


The  solution  can  now  be  completed  in  precisely  the  same  manner 
as  before.  This  method— of  examining  the  detailed  behavior,  possibly 
to  several  terms,  in  the  expansion  of  Y^(s)  as  s in  R^  —is 
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frequently  used  to  determine  unknown  constants  arising  from  boundary 
conditions.  Another  common  method,  which  Is  Illustrated  In  the  next 
section  on  waves  In  beams,  Involves  obtaining  a solution  like  (5.6) 
for  general  values  of  the  unknown  constants,  and  then  arguing  that 
unless  the  constants  have  certain  special  values,  a plus  function  will 
have  a pole  somewhere  in  or  a minus  function  will  have  a pole  some- 
where in  R . In  complicated  problems  involving  coupled  elastic  plate/ 
acoustic  fluid  motions  the  plus  or  minus  function  which  apparently 
has  the  singularity  may  not  be  the  obvious  or  Y function,  but  some 
more  complicated  related  function. 
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6.  WAVES  ON  BEAMS 


We  now  replace  Che  strings  of  ^2  by  beams  of  specific  mass 
and  bending  stiffness  B^,  B^  in  x < 0,  x > 0,  respectively.  The  free 
wavenumbers  will  be  denoted  by  as  before,  for  time  dependence 

exp(-ia)t),  where 

kj  - (py/B^)  , k;;  - (PjwVBj).  (6.1) 


For  the  moment  we  leave  conditions  at  the  junction  x “ 0 unspecified 

and  Cake  a wave  with  displacement  exp(lk  x)  incident  from  x > -<*>, 

o 

denoting  the  total  displacements  again  by  y + exp(ik  x)  in  x < 0 
o 

and  by  y in  X > 0.  This  ensures  exponential  decay  of  y(x)  as 

X -►  + oo  and  as  X “ provided  Im  k , Im  k are  given  small  positive 

o 1 

values,  fcr  we  anticipate  that  in  x > 0 


y(x)  » A exp(lkjX)  + B exp(-kjX) 

and  (6.2) 

y(x)  “ C exp(-ik  x)  + D exp(k  x) . 

o o 

The  near-field  terms  here,  B exp(-kjX)  and  D exp(k^x),  decay 

as  X -►  + oo  and  as  x -►  - <»,  respectively,  because  k and  k,  have  positive 

o I 

real  parts;  they  arise  because  the  equations  of  motion  are 


^ -k'*y“0  inx<0 

dx-  o 


k-  y 
dx-  • ^ 


0 


in  X > 0 


Taking  half-range  Fourier  transforms  of  (6.3)  gives 


(6.3) 
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(s'*  - ki)  Y^(8)  - y'"((H)  - isy"  (CH-)  - s^y'COf)  + i8*y((H) 


(6.4) 


(s'*  - k'*)  Y (s)  - -y"'(0-)  + i8y"(0-)  + s^y*  (0-)-ls®y (0-)  (6.5) 

o — 

the  first  of  these  holding  in  R^,  the  second  in  R . Now  the  poles 
s “ +kj,  s ■ +ikj  lie  in  R^  (see  Figure  3),  so  that  Y^(s)  will  have 
poles  at  +kj,  + ik^  unless  we  choose 

y"’(Of)  - ikj  y"((>+)  - kj  y’(0+)  + ikj  y((H)  « 0 (6.6) 

y"’((H)  + k,  y"(CH-)  + kj  y'(0+)  + kj  y(0+)  » 0 (6.7) 

Similarly,  Y (s)  will  have  poles  at  s » ® ^ unless 

-y"'(0-)  - Ik  y"(0-)  + k^  y'(O-)  + ikV(O-)  - 0 (6.8) 

o o o 

-y"'(0-)  + k y"(0-)  ^ k^  y’(O-)  + k^  y(O-)  = 0 (6.9) 

o o o 

Whatever  the  boundary  conditions  at  x = 0,  conditions  (6.6-6. 9)  must  be 
satisfied.  (Two  analogous  relations  could  have  been  deduced  in  ^2,  but 
it  seemed  unecessary  to  emphasize  such  a point  at  an  early  stage.) 

Four  further  conditions  may  be  Imposed  at  x * 0,  at  least  one  of 
these  being  a nonhomogenous  condition,  so  that  we  shall  have  a uniquely 
soluble  set  of  eight  equations  for  the  eight  unknown  boundary  constants. 

We  shall  not  take  any  particular  set  of  boundary  conditions,  as  these 
tend  only  to  lead  to  complicated  expressions  without  any  special  structure. 

By  addition  of  (6.4)  and  (6.5)  we  get 

(s'*  - k*;)  Y_^(8)  + (s"  - k^)  Y_(8)  - Qj(8)  + Q^(8)  (6.10) 
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where 


Qj(s)  ■ is*y(0+)  - 8*y' (Of)  - i8y"(0<-)  + y'"(0+) 

Q^(s)  --isV(O-)  + s*y'(0-)  + i8y"(0-)  - y’"(0-) 

and  conditions  (6. 6-6. 9)  are  aatlefied,  80  that 

Qj(kj)  - Q/ik,)  - 0 
o o o o 

Equation  (6.10)  holds  in  the  atrip 

D:  -k, .<  Im  s < + k . 

11  oi 

For  simplicity  we  can  take  k^^  = kj  so  that  the  strip  D is  symmetric 

about  the  real  axis  and  the  points  +k  , -He  lie  on  the  upper  boundary 

o * 

of  D,  the  points  ~k^,  -k^  on  its  lower  boundary.  Because  k^^  is 
supposed  to  be  very  small,  the  other  points  of  Interest,  ’♦‘ik^,  -Hkj, 
and  ~i-kjj»  ~ikj,  lie  well  above  and  well  below  D,  respectively.  If 
therefore  we  work  in  the  interior  of  D we  can  divide  (6.10)  through 
by  (s'*  - k^)  say,  and  get  a W-H  equation 

Qi(s)  + Qo(s) 

K(s)  Y^(s)+Y_(8)  - — — (6.13) 

(s  - k^) 

in  which  the  kernel  is 

s'*  - k** 

K(8)  - — J (6.14) 

® ” *^0 

The  U-H  product  split,  into  factors  analytic,  non-zero  and  of 
algebraic  behavior  at  infinity  in  , respectively,  is  again  obvious: 


(6.12) 


(6.11) 
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and  we  have  to  make  an  additive  split  of  the  right  hand  side.  Now 
because  of  (6.12),  Qq(s)  must  contain  the  factors  (s  + » 

while  Qj(s)  must  contain  the  factors  (s  - k^Xs  - ik^)  so  that  we 
can  write 


Q (s)  = (s  + k ) (s  + ik  ) (a  s + b ) 
Q o o o o 

Q^(s)  E (s  - k^)(s  - lk^)(a^s  + b^) 


(6.17) 


where  the  coefficients  a , a , b , b.  are  known  when  any  particular 

o 1 o 1 

set  of  conditions  Is  specified  at  x ~ 0.  Now  the  right  side  of  (6.16) 
has  the  form 

as+b  . a,s+b, 

o o + t 1 

(s  - kj)(s  - ikj)  (s  + k^)(s  + ik^) 


which  is  already  in  the  desired  form  G (s)  + G^(8).  The  reason  for 
this  Is  that  analytlclty  arguments  have  already  been  used  to  remove 
pole  singularities  where  they  are  not  permitted  and  because  pole 
singularities  are  the  only  kinds  of  singularity  which  are  present  In 
these  one-dlmenslonal  problems  this  effectively  means  that  the  additive 
split  must  already  have  been  carried  out. 
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By  the  usual  arguments  we  then  have  the  solutions 


K^(s)  Y^(s)  - G^(s)  - P(s) 


- Hi}  ■ 


(6.18) 


where  P(s)  is  a polynomial.  And  because  G^(s)  * 0(s“‘),  K^(s)  - 0(1) 
and  Y^(s)  - 0(s”')  (because  y is  finite  at  x = 0)  at  Infinity  in  their 
respective  half-planes  of  analyticity,  Che  polynomial  P(s)  must  be 
Identically  zero,  so  that 


a^s  + bi 


Y r ® 

+'  (s  + kj ) (s  + Ikj ) 


Y (s) 


a s + b 
o o 


(s  - k ) (s  - ik  ) 
o o 


The  inverse  transforms  can  then  be  performed  explicitly,  the  poles 

s “ k and  s « ik  in  R,  gives  rise  to  the  reflected  wave  and  a 
o o + 

decaying  mode  in  x 0,  the  poles  s = -kj  and  -ikj  gives  the  trans- 
mitted wave  and  a decaying  mode  in  x > 0.  Specifically  we  find 


(a  k,  - b , ) 

y(x  > 0)  = i ^ ~ exp(ikjx) 


(ajiki  - bj) 

^ (i"  i) — exp(-kix) 


(a  k + b ) 

y(x  < 0)  » 1 — — exp(-lk^x) 


k (1  - 1) 
o 


(a  ik  - b ) 

^ “kV”  i-f- 


which  is  of  the  anticipated  form  (6.2).  For  any  specified  boundary 
conditions  Che  values  of  y and  its  first  three  derivatives  are  known 
at  X ■ 0^,  so  that  the  coefficients  a^,  b , a^,  b^  in  (6.20)  can  be 


(6.19) 


(6.20) 


: I 
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I 


obtained  from  a direct  comparison  of  (6.17)  with  (6.11).  But  of 
course,  the  W-H  was  never  intended  to  be  used  for  solving  such  simple 

! 

problems  as  these  one-dimensional  ones.  We  are  using  them  as  the 
simplest  vehicle  on  which  many  features  of  common  occurrence  in  the  W-H 
I method  can  be  demonstrated. 
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7.  GENERALIZED  W-H  EQUATIONS:  THREE-PART  BOUNDARY  VALUE  PROBLEMS 


The  standard  W-H  equation,  (2.26),  arises  In  many,  though  by  no 
means  all,  boundary  value  problems  In  which  different  boundary  data  are 
prescribed  on,  say,  x < 0 and  x > 0.  Many  problems  of  Interest  In 
acoustics  fall  Into  this  category:  for  Instance,  the  problem  of  energy 
conversion  from  the  elastic  to  the  acoustic  mode  when  a surface  wave 
In  an  elastic  plate  encounters  a junction  In  the  plate  across  which 
the  plate  properties  change  abruptly  can  be  modeled  In  terms  of  two 
semi- Inf Inlte  plates  x < 0 and  x > 0 for  many  purposes.  The  boundaries 
concerned  do  not  always  have  to  coincide  with  just  the  x-axls.  For 
example,  the  standard  W-H  equation  arises  In  the  diffraction  of  waves 
by  an  open-ended  parallel  plate  waveguide,  or  by  an  open  ended  circular 
duct,  provided  these  are  both  seml-lnf Inlte.  On  the  other  hand,  dif- 
fraction by  three  parallel  equl-spaced  seml-lnf inlte  plates  Is  a com- 
pletely open  problem,  though  diffraction  by  an  Infinite  cascade  of 
seml-lnf inlte  staggered  plates  is  a relatively  simple  standard  W-H 
problem.  If  the  waveguides  referred  to  above  have  closed  ends  (diffrac- 
tion by  a seml-lnf Inlte  thick  rigid  plate,  or  by  a seml-lnf Inlte  solid 

i 

I 

rod)  the  standard  W-H  method  does  not  lead  to  a closed  form  solution, 
but  to  an  Infinite  set  of  coupled  linear  equations,  whose  solution  can 
only  be  approximated  by  the  solution  of  a finite  subset  of  the  equations 
In  the  low  frequency  limit. 

Thus  It  Is  clearly  difficult  to  give  any  general  guidelines  as  to 
when  the  standard  W-H  method  would  work  except  to  say  that  It  will  not 
work,  without  modification.  In  the  case  of  three-part  boundary  value 
problems  where  data  Is  given  on,  say,  - «><  x<0;  0<x<il; 


32 


£ < X < + « in  three  different  forms.  There  are,  however,  a number  of 
problems  which  have  been  successfully  Cackled — though  always  only  in 
approximate  form — in  this  category,  using  modified  W-H  methods.  We  quote 
as  examples  the  problems  of  diffraction  by  an  infinite  rigid  strip  of 
finite  width  in  the  high  frequency  limit,  of  the  resonances  of  a circular 
tube  of  finite  length  open  at  both  ends  and  of  wave  motion  in  an  elastic 
plate  set  in  an  infinite  rigid  surrounding  baffle. 

Now  in  many  of  these  problems  Che  approximate  solution  Is  achieved 
on  Che  basis  of  a "weak  Interaction"  simplification.  For  example.  In 
diffraction  of  waves  of  wavenumber  k^  by  a rigid  strip  of  width  £ in 
the  high  frequency  limit,  >>  1*  ofie  can  suppose  that  to  first 

order  each  edge  Is  unaware  of  the  presence  of  the  ocher,  so  that  one 
can  start  the  approximation  with  two  semi-infinite  problems  of  standard 
W-H  type.  Then  if  the  incident  wave  amplitude  is  0(1),  the  first  inter- 
action of  one  edge  with  the  other  will  be  through  a cylindrical  wave 

emanating  from  one  edge,  due  to  the  incident  wave,  and  of  amplitude 
-± 

0(k  £)  near  the  other  edge.  For  the  second  approximation  we  there- 
o 

fore  solve  another  two  semi- inf inite  problems,  but  now  with  more 
complicated  forcing  terms  arising  from  the  mutual  interaction  between 
the  edges.  While  one  can  see  how  to  continue  the  process  in  that 
simple  case,  it  Is  advantageous  to  derive  a modified  W-H  equation  whose 
approximate  solution  throws  up  these  successive  interaction  problems  in 
a natural  way.  The  advantage  is  that  one  can  see  from  the  generalized 
W-H  equation  what  to  try  in  more  complicated  problems  where  the 
physical  situation  is  less  clear.  For  example,  near  resonance  the 
"weak  interaction"  sort  of  approximation  is  quite  Inappropriate,  as  the 
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whole  phenomenon  depends  on  strong  coupling  between  the  ends  of  the 
system. 

We  shall  try  here  to  use  our  one-dimensional  examples  to  illustrate 
the  possibility  of  generalizing  the  W-H  technique  to  deal  with  three- 
part  problems.  For  such  one-dimensional  problems,  however,  there  is 
not  generally  any  "weak  interaction"  approximation  that  one  can  make, 
because  waves  on  strings  do  not  decay  in  the  way  that  two  and  three 
dimensional  acoustic  fields  do,  so  that  in  a sense  one  is  always  con- 
fronted with  the  strong  coupling  situation.  Nonetheless,  a number 
of  interesting  points  arise  in  the  string  problems  which  have  direct 
comparisons  in  more  serious  three-part  problems. 

Consider  a uniform  string  of  line  density  in  - <»  < x < 0 
and  in  il  < X < + “,  these  two  semi-infinite  strings  being  joined  by  a 
string  of  line  density  in  0 < x < i.  A wave  with  displacement 
exp  ik^x  is  Incident  from  x = - «>,  and  we  want  to  find  the  displacements 
everywhere,  subject  to  conditions  at  x = 0 and  at  x = £ which  we  will 
specify  later. 

Write  the  total  displacement  in  x < 0 as  y + exp  Ik  x,  so  that  y 

o 

represents  an  outgoing  field  as  x -►-<».  Define 

Y_(s)  = J y(x)  exp  isx  dx  (7.1) 

so  that  Y (s)  is  analytic  in  R (Im  s < k .)  and  of  algebiaic  decrease 
- - ol 

everywhere  at  infinity  in  R (because  y(x)  is  finite  at  x “ 0) . From 
the  differential  equation 

+ k*y  ”0  (-  <»  < X < 0) 

dx*  o 


34 


we  have 


(s^  - k^)  Y (s)  = + y'(O-)  - isy(O-) 
o - 


and  in  order  to  avoid  having  a pole  at  s “ -k^eR_we  must  have 


+ y’ (0-)  + ik^  y(O-)  = 0 


(7.2) 


(7.3) 


Also  as  a check,  as  | s | “ in  R , 


Y_(s)  ~ y(O-) 


which  we  also  get  from 
Y (s) 


-r 


^ . N isx  . 

y (0-)  + . . . ) e dx 


using  Watson's  lemma. 

The  same  differential  equation  holds  in  < x < <».  Define 


'+<=)  ■ r 


) exp  isx  dx 


(7.4) 


(7.5) 


y(x)  being  the  total  displacement  in  £ < x < Then  Y_^(s)  is  analytic 
in  the  upper  half-plane 


R : Im  s > - k . 
+ oi 


because  of  the  exponential  decrease,  exp(-k^^x)  of  y(x)  as  x ■►  + “. 
However,  Y_^(s)  does  not  now  have  algebraic  behavior  at  infinity  in  R_^. 
To  find  the  behavior  we  write 
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Y+(s)  = r y(X  +1)  dX 

''  o 

= y(X  +i)  dX 

• /^ 


~ e 


ls£. 


f 

• fy 


[y(Jl+)  + X y'(U)  + 


, IsX  y 
. . j e dX 


^ i Isa  ,0^,  ^ 
= + - e y(a+)  + 


•(■?) 


as  s -*•  “ in  R^.  Thus  Y^(s)  has  exponential  behavior,  exp  is^,  at 
infinity  in  R^.  From  the  differential  equation  we  have 

(s'  - k')  Y (s)  = [-y'(a+)  + is  y(a+)]  e^®^ 
o + 


(7.6) 


(7.7) 


from  which  (7.6)  is  obvious.  Again,  there  is  a pole  at  s ” k eR^  unless 

o + 


-y'(a+)  + ik^  y(a+)  = o 


For  the  middle  portion  of  the  string  we  define 

Yj(s)  = I y(x)  exp  isx  dx 
''0 

Since  the  integration  is  over  a finite  range  only  and  the  integrand  is 
bounded  and  continuous,  Yj(s)  is  analytic  in  the  entire  complex  plane 
(the  Fourier  Transform  of  a function  with  compact  support  — i.e., 
vanishing  outside  some  finite  range,  is  an  ENTIRE  function).  As  s « 

with  Im  s > 0, 


(7.8) 


(7.9) 


Yj(8  ) ~ r (y((H-)  + xy’(0+)  + ...)  e""®*  dx  - 0(s~‘) 
•'0 


isx 


(7.10) 


because  when  Im  s > 0 only  small  values  of  x contribute  to  the  Integral, 
so  that  we  can  expand  y(x)  about  x ^ 0+  and  also  extend  the  integration 


36 


up  to  As  s » with  Im  8 < 0 we  need  to  write  x = £ - z,  to  get 

Yj(s)  = exp(is£)  I y(£  - z)  exp(-isz)  dz 
•'o 

and  because  of  the  exponential  exp(-isz)  we  can  expand  about  z > 0 again 
to  get 


Yj(s)  - exp  is£  j [y(£-)  - zy'(£-)  + ...]e 

• r\ 


dz 


= exp(is£)  — y(£-)  + 0(-^) 
s s 


(7.11) 


Thus  the  entire  function  Yj(s)  is  algebraically  small,  0(s“M  or  smaller, 
in  Im  s > 0,  but  exponentially  large  like  exp(isl!,),  in  Im  s < 0. 

The  differential  equation 


^ + k^y  = 0 
dx2  ' 

(where  y is  the  total  displacement  in  0 < x < £)  gives 
(s^  - k^)  Yj(s)  = (y*(£-)  e^®^  - y'(0+)] 

- is[y(£-)  e^®^  - y(0+)]  (7.12) 

and  confirms  (7.10)  and  7.11).  Further,  Yj(s)  can  have  no  singularities 
for  any  finite  value  of  s,  so  that  the  right  side  of  (7.12)  must  vanish 
for  both  s “ k^  and  s »•  -k^ , giving 

(y'(£-)  - y'((H-)]  - ikj[y(il-)  - y(0+)]-  0 , (7.13) 

ly'(£-)  - y'((H)]+  ikj(y(£-)  - y(0+)]  - 0 . (7.14) 
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To  make  the  algebra  minimal  we  choose  a definite  set  of  conditions 
at  X ~ 0 and  at  x ■ namely  the  simple  junction  conditions  that  there 
is  no  change  to  the  total  displacement  or  to  the  slope  at  x ■ 0,  x * 1, 
the  tension  in  all  three  strings  being  the  same.  Thus  we  take 


y(O-)  + 1 = y(0+) 

y(£-)  = y()l+) 

y'(O-)  + Ik  - y'((H-) 
o 

y'()l-)  = y'(£+) 


(7.15) 


If  k k,  the  set  of  8 equations  (7.3,  7.8,  7.13,  7.14,  7.15)  has  a 
o » 

finite  solution  provided 

sin  kj£  ^ 0 (7.16) 

(excluding  resonance  of  the  middle  portion)  and  then 


yW 


y(CH-) 


Zlkoki 

k^-k? 
o 1 


co9ec(kj£) 


2k^ 

o 

k^-k^ 
o 1 


2ik  k, 

— cot  k £ 

k^-k"  ' 

o 1 


(7.17) 


etc.  This  of  course  completes  the  solution  for  these  simple  problems, 
for  now  that  all  constants  are  known,  T^(s)  is  known  from  (7.7), 


Y_(8)  from  (7.2)  and  Yj(s)  from  (7.12). 

A method  which  shows  how  a generalized  W-H  equation  may  be  treated 
ignores  the  detailed  solution  (7.17),  and  instead  eliminates  the  unknown 
constants  from  (7.7),  (7.2),  and  (7.12)  in  just  the  same  way  that  the 
corresponding  constants  were  eliminated  in  ^2  to  get  a standard  W-H 
equation.  Thus  we  write 


38 


(s*  - kj)  Yj(s)  - e^®*'[y'(il-)  - Is  y(£-)] 
- ly’((H-)  - is  y(CH-)] 


(7.18) 

(7.19) 

This  is  the  required  generalized  W-H  equation,  a single  equation  for 
three  unknown  functions  Y^(s),  Y (s)  and  Yj(s),  given  the  kernel  K(8) 
and  the  forcing  field  (s  + k^)~*. 

Equations  of  this  kind  have  been  considered  by  Noble  [5  p.  196]. 

Methods  exist  for  solving  such  equations  approximately  in  "high 
frequency"  limits  in  which  the  finite  part  of  the  boundary  is  many 
wavelengths  long,  in  some  appropriate  sense.  Generally  these  methods 
rely  on  weak  interaction  between  the  ends  x ■ 0 and  x > though  as 
remarked  before,  the  ends  are  rather  delicately  coupled  in  cases  such 
as  the  resonance  of  a finite  open-ended  waveguide,  and  methods  have 
also  been  developed  to  deal  with  such  cases.  Here,  because  only  pole 
singularities  are  involved,  it  is  possible  to  solve  (7.18)  exactly. 


- (s^  - k*)  Y (s) 
o + 


- (s^  - kM  Y (s) 
o - 


+ i(s  - k ) 
o 


on  use  of  (7.15),  (7.2),  and  (7.7),  so  that 


Y^(s)  + Y_(s)  + K(s)  Y^(s) 


8 + k 


where  the  kernel  is,  as  in  02, 

s^  - k? 

K(s)  = 


8^  - k^ 
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We  carry  the  analyses  through  In  a general  form  as  far  as  possible,  to 
Indicate  the  procedure  which  has  to  be  followed  In  more  complicated 
problems. 

First  of  all  we  extract  a factor  exp  isli  from  Y^(8) , writing 


Y^(s)  - exp(is«,)  Z^(s) 


(7.20) 


so  that  Z^(s)  » 0(s  *)  at  infinity  in  R^.  The  necessity  for  doing 
this  will  be  apparent  in  a moment.  Then  write  K(s)  “ K^(s)  K (s)  as 
usual,  and  divide  through  by  K (s)  to  get 


ls£_  , . 
e Z^(s) 

' K (s) 


Y_(s) 


i 

(s  + k ) K (s) 
o — 


The  second  term  on  the  left  is  analytic  in  R and  0(s“*)  at  infinity 
there,  the  third  is  analytic  in  R^  and  0(s~*)  at  infinity  there.  The 
term  on  the  right  can  be  split  in  the  familiar  way  as 


i 

(s  + k ) K (s) 
o - 


G^(s)  + G_(s) 


(7.21) 


and  we  also  make  an  additive  split  of  the  first  term  in  this  way, 
isJl 


Z^(s) 


K (s) 


= U^(s)  + U_(s) 


(7.22) 


using  a general  theorem  to  be  given  in  a moment.  G^(s)  are  each 
0(8~^)  at  infinity  in  R^,  and  we  assume  that  is  also  true  of  U^(s). 
Then  we  have 


U^(8)  + K^(8)  Y (s)  - G^(8) 


- G (8) 


Y-(s) 

K_(8) 


- U_(8) 
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and  each  side  Is  the  representation.  In  or  R as  the  case  may  be, 
of  a single  entire  function  which  behaves  like  s~*  everywhere  at 
Infinity,  and  Is  therefore  Identically  zero.  Hence 


K^(8)  Y^(s)  - G^(8)  + U^(8)  - 0 
Y_(8) 

- FTI)  - ° • 


(7.23a) 


(7.23b) 


Now  return  to  the  generalized  W-H  equation,  and  this  time,  instead 
of  dividing  by  K (s),  we  divide  by  exp(is£,)  K^(8)  to  get 


Z^(s)  XJs) 

K^(s)  K^(s)  ® 


K_(s)  Yj(s) 


ie 


-isH 


(s  + k ) K,  (s) 
o + 


(7.24) 


The  first  term  is  analytic  in  R and  0(s"*)  at  Infinity  there  while 
the  third  is  analytic  in  R and  0(s“')  there  only  because  of  the  factor 
exp(-is£.).  If  we  had  not  divided  by  exp  isil  we  would  have  been  left 
with  a third  term  which  was  analytic  in  R , but  exponentially  large 
at  infinity  in  R (see  Eq.  7.11)  and  the  function  theoretic  argument 
would  not  go  through  (in  particular,  the  entire  function  would  not  be 
zero,  but  would  be  exponentially  large  in  R , and  there  is  no  general 
way  of  constructing  functions  of  this  kind) . 

The  function  appearing  second  on  the  left  in  (7.24)  has  mixed 
properties,  so  we  again  try  to  split  it  as 


-isJl 

e 


Y (8) 


K^(s) 


- V^(8)  + V_(8) 


(7.25) 
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with  v^(8)  - 0(8“ ‘)  at  infinity  in  R^.  On  the  other  hand,  the  function 
occurring  on  the  right  of  (7.24)  is  a@functlon — but  it  i8  exponentially 
large,  like  exp  (Sjf.)  in  the  upper  half-plane,  and  80  it  i8  neceasary 
to  make  a split 


. -is^ 

(s  . k ) K,(s)  ■ ^ "-<*> 


(7.26) 


with  H^(s)  • 0(s“*)  in  R,.  in  order  to  remove  this  exponential  increase. 

Now  we  can  split  (7.24)  into(5 andQ parts , each  of  which  is  0(s  *) 
at  infinity  in  R+  and  each  of  which  therefore  vanishes  identically. 

Thus  we  get 


H_(s)  - V_(s)  - e“^®^  K_(s)  Yj(s)  - 0 . 


(7.27a) 


Z.(s) 

•n-T-r  + V (s)  - H (s)  - 0 . (7.27b) 

K^(s)  + + 

The  equations  (7.23a,b;  7,27a,b),  obtained  by  the  W-H  argument,  are  not, 
in  general,  solutions  to  the  problem,  but  they  constitute  a pair  of 
integral  equations  which  are  in  a form  suitable  for  approximate  solution 
in  the  "high  frequency  limit"  (l.e.,  for  k^*.  » 1,  kjf  » 1 in  general). 

To  see  this  we  have  to  use  the  general  formula  (see  Section  10)  for 
expressing  a function  F(8),  analytic  in  the  strip  D and  with  suitable 
behavior  at  infinity  in  D,  as  the  sum  of  functions  analytic  in  R^  and 
0(8  ')  at  Infinity  there.  If 

F(8)  - F^(8)  + F_(s) 

where  the  functions  are  to  have  the  stated  properties,  then 
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dt 


''7.28a) 


F (3) . ^ f .mi 

^ 2Tri  J t - 


F (8) 


zL.  f m)  dt 

2Trl  J t - 8 “ ’ 


(7 


where  in  (7.28a)  the  path  runs  from  - to  + “ in  the  strip  D passing 
below  the  point  t = s while  in  (7.28b)  the  path  passes  above  t * s. 

Applying  (7.28)  here  we  have,  for  (7.23b)  and  (7.27b)  in  particular, 
Y (s) 


Z^(t) 


2iTi  J K^(t)(t  - s)  ° ^ 


Z (s)  1 f e 

IT 


Y (t) 


(t)(t  - 8) 


dt  « 0 , 


(7 


(7 


the  forcing  fields  G_(s)  and  H^(s)  being  known,  in  principle.  Clearly 
(7.29),  (7.30)  are  a pair  of  coupled  Integral  equations  for  the  unknown 
functions  Y_(s)  and  Z^(s)  which  determine  the  reflected  and  transmitted 
fields  in  X < 0,  X > 1,  respectively.  Once  these  are  known,  the  field 
Y^(8)  in  the  middle  portion  of  the  string  can  be  found  directly  from 
the  generalized  W-H  equation  (7.18),  for  example. 

To  see  the  structure  of  these  equations,  suppose  that  the  Integral 
term  in  (7.29)  were  zero.  Then  we  would  have 


Y_(8)  - K_(8)  G_(8) 

and  noting  the  definition  (7.21)  of  G (s)  and  comparing  with  (3.11) 
we  see  that  this  Y_(8)  is  precisely  the  field  in  x < 0 if  the  wave 
exp(ik^x)  were  Incident  upon  a semi- infinite,  rather  than  finite, 
string  to  the  right.  Similarly,  the  situation  Z^Cs)  ■ K^(8)  H^(8) 
is  the  solution  for  a semi-infinite  problem  of  reflection  at  the 


.28b) 


.29) 


.30) 
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junction  x ■ the  left  hand  portion  of  the  string  now  being  regarded 

as  extending  to  x ■ - The  Integral  terms  represent  Interactions 

between  the  ends  x ■ 0,  x •*  H,  giving  rise  to  a sequence  of  reflections 
and  transmissions.  Methods  have  been  devised  for  dealing  approximately 
with  these  Interactions  both  when  they  are  weak  and  when  they  are  strong, 
as  In  resonance  situations  (see  Noble  [5  [Ch.5]). 

Here  It  Is  of  course  possible  to  solve  the  coupled  integral  equa- 
tions completely.  First,  however,  we  show  how  they  may  be  uncoupled 
when  the  kernel  K(s)  Is  even  in  s,  as  It  Is  here.  It  Is  then  possible 
to  work  In  a strip  D which  is  completely  symmetrical  about  the  real  axis, 
so  that  Y (-s)  is  a (^function  with  the  same  domain  of  analytlclty 
as  Z^(s),  while  Z^(-s)  is  aQ function.  It  also  follows  from  Sec- 
tion 10  that  K^(-s)  ■ K_(s)  for  sCD. 

Consider  then  (7.29)  and  (7.30)  fcr  seD,  and  choose  the  path 
— — above  s to  be  a straight  line  from  +ia  to  +“  +ia,  a > 0. 
Choose  — to  be  the  Image  of  this  path  in  the  real  axis,  going 
from  -la  to  +«-ia.  In  (7.29)  change  s to  -s;  the  integration  path 

can  be  chosen  to  run  above  both  s and  -s,  and  so 

^Itl 

+ I 5 

,,  y-iTi  I V 

+ 


Y_(-s)  1 f e^ 

■ K (s)  2Tri  J F 


Z+(t) 
(t) (t+s) 


dt 


(7.31) 


In  (7.30)  change  the  Integration  variable  from  t to  -t;  (-t)  runs 
from  +ia  +«oto  +ia  -«>as  t runs  from  -la  -<»  to  -la  +<*>  , so  that  the 
path  for  (-t)  is  the  same  as  that  In  (7.31).  This  gives 

.ltlY_(_t)  dt 


. Ms)  1 f e"'-*Y_(-t) 
-H^(8)  + - 2^^  J K^(t)(t+n) 


- 0 , 


(7.32) 
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Now  add  (7.31)  and  (7.32)  defining  the  difference 
D^(8)  - Z^(8)  - Y_(-8)  , 

to  get 

iti. 

D (8)  if® 

[G_(-8)  -H^(s)]  + K^(s)  ^ J K_(t)(t+8) 

while  if  we  eubtract  (7.31)  from  (7.32)  and  define 


S^(8)  - Z^(8)  + Y_(-s) 

we  get 

S (s)  /.  S^(t)  dt 

(-G_(-s)  -H_j^(s)]  K^(s)  " 2kI  J K_(c)(t  + s)  ° 

80  that  (7.33)  and  (7.34)  are  a pair  of  similar  uncoupled  integral 
equations. 

Let  us  now  look  at  the  forcing  functions  in  these  equations. 
The  additive  split  of  G(s)  in  (7.21)  is  simple,  giving 


G^(s) 


i 

(s  + k ) K (-k  ) 
0 — 0 


G (s) 


i 

8 + k 

o 


1 

K (s) 


\ 


> 


/ 


For  the  function  H(8)  in  (7.26)  we  have  to  uae  the  Cauchy  Integrala, 
which  give 


H^(8) 


2iii  J (t  + k ) K.  (t)  Ct  - s) 

• o 


(7.33) 


(7.34) 


(7.35) 


(7.36) 


We  cannot  complete  the  contour  with  a large  semicircle  in  R^,  because  the 
factor  exp (-iti)  is  exponentially  large  there.  Instead  we  complete  the 


contour  with  a large  semicircle  In  R_,  along  which  the  Integral  vanishes 
because  of  exponential  smallness  of  exp(-lt£)  there.  To  Interpret  the 
meaning  of  K^(t)  In  R_  we  write  It  as  K(t)/K  (t),  so  that 


u , , 1 f le~^*^^  K_(t) 

”+^®^“27rl  J (t  + k ) K(t)(t  - s) 

V o 

The  pole  at  t = s lies  outside  the  contour,  K (t)  Is  analytic  within 
the  contour,  and 

(t  + k^)  K(t)  = (t  + kj) 

so  that  the  only  singularity  within  the  closed  contour  Is  at  t = -k^ , 
and  therefore 


(7.37) 


H^(a) 


2tt1 


(-2Trl) 


4 If  0 

le^*^!  K_{-kj) 


, lk,«. 
le  1 

(s  + kj)  ’ 


(7.38) 


H (s)  - 


le 


-isi 


(s  + k^)  K^(s) 


, ik,«. 
le  1 

(s  + kj) 


8 + k 


-(e 


-IsH 


Ik.l!., 

- e 1 ) 


(7.39) 


This  Instructive  example  shows  the  Importance  of  removing 
exponential  growth  at  Infinity;  although  the  original  H(s)  Is  analytic 
In  R^  It  Is  not  algebraically  small  at  Infinity  there,  and  use  of  the 
Cauchy  Integrals  shows  how  It  can  be  split  Into  H^(8)  + H (s) , with 
H^(8)  - 0(s"*)  at  Infinity  In  R^  as  In  (7.38)  and  (7.39). 

Coming  now  to  the  Integral  terms  In  the  Integral  equations  (7.33), 
(7.34),  we  can  this  time  complete  the  contour  with  a large  semicircle 
In  R^,  because  along  that  semicircle  the  factor  exp(ltil)  will  be 
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exponentially  small  and  Che  Integral  along  the  semicircle  will  contribute 


(7. AO) 


(7.41) 


nothing.  The  functions  D^(t),  S^(t)  are  analytic  within  the  closed 

contour,  and  the  function  1/K_(t)  ■ K^(t)/K(t)  has  a pole  at  t “ only. 

Therefore  (7.33)  gives 

D.(s)  (k,  - k ) D.(k.) 

[G.(-s)  -H^(s)]  0 

and  (7.34)  gives 

S.(s)  (kj  - k ) S.(k.) 

[-G_(-s)  .H^(s)]  + ^ i—  . 0 

The  unknown  constants  D_j^(kj),  S^(kj)  are  found  by  putting  s * k,  in 
(7.40),  (7.41)  and  then  we  have  completely  determined  the  functions 
D_^(8),  S^(s),  from  which  Z^(s) , Y (-s)  can  be  found,  and  hence  the 
whole  solution  is  determined.  While  there  are  no  difficulties  of 
principle,  the  algebra  Is  very  tedious  and  there  Is  no  point  in  giving 
It  here. 

The  only  remaining  point  of  Interest  concerns  the  Inversion  of 
the  Fourier  Integral  for  this  three-part  problem.  We  have 

y(x)  - j [Y_(8)  + Yj(8)  + Y^(s)]  exp(-i8x)  ds  (7.42) 

with  the  Integral  along  a path  from  -<»co  -l-<»ln  D,  and  we  recall  that 


Y (s)  - 0(8-*) 


Yj(8) 


Y^(8) 


0(8-*) 
Isl 

-0(^) 

Isl 

- 0(^) 


as  8 -►  0“  in  R 


as  8 -►  » In  R. 


as  8 00  in  R 


as  8 -►  00  in  R. 


\ (7.43) 
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that  Y^(8)  are  analytic  in  R^,  respectively,  and  that  Y^(8)  is  analytic 
everywhere  (an  entire  function).  For  x < 0,  deform  the  contour  into  the 
upper  half-plane  R^.  The  contributions  from  the  large  semicircle  in  R^  to 
the  Yj  and  Y^  Integrals  are  zero,  because  exp(-lsx)  is  exponentially  small 
and  Yj  and  Y^  are  at  least  as  small  as  s ' at  infinity  in  R^.  Further, 

Y^  and  Y^  are  analytic  within  the  closed  contour  , and  so  their 


contour  integrals  vanish.  Thus  for  x < 0 


(s)  exp(-isx)  ds 


(7.44) 


which  we  may  try  to  evaluate  by  completing  the  contour  in  R^ — thougn 
note  that  we  know  nothing  in  advance  about  the  behavior  of  Y (s)  at 
infinity  in  R^. 

If  X > £.  we  complete  the  contour  with  a large  semicircle  in  R , 
along  which,  although  Y^(s)  is  exponentially  large,  the  product 
Y^(s)  exp(-isx)  is  exponentially  small  when  x > £.  Similar  arguments 


then  give 


- f Y 
2tt  J + 


(s)  exp(-isx)  ds 


(7.45) 


Finally,  if  0 < x < £,  close  the  contour  for  Y (s)  in  R , that 
for  Y^(s)  in  R_j^,  and  we  find  that  the  integrals  of  Y (s) , Y^(s) 


vanish,  leaving 


2v  J 1 


(s)  exp(-isx)  ds 


(7.46) 


Because  of  (7.43)  it  is  only  possible  to  close  the  contour  here  in  R , 
Yj(s)  exp(-isx)  being  exponentially  large  in  R^  when  0 < x < £. 

These  arguments  apply  generally  in  three-part  boundary  value 
problems.  Here  the  integrals  (7.44-7.46)  can  all  be  evaluated  by 


J 


residue  calculus,  since  only  simple  poles  are  Involved,  leading  again 
to  results  which  can  be  confirmed  using  elementary  methods. 


The  purpose  of  this  section  has  not  been  to  solve  a particular 
string  problem,  but  to  Illustrate  how  the  W-H  method,  applied  to  three- 
part  problems,  leads  In  general  not  to  a solution,  but  to  a pair  of 
coupled  Integral  equations.  For  an  even  kernel  K(s)  we  can  decouple  the 
equations,  and  deal  with  a pair  of  similar  Independent  Integral  equations 
for  the  functions  D^(s) , S^(s) . In  our  case  the  Integral  terms  can 
actually  be  evaluated  explicitly  in  terms  of  unknown  constants  D^Ck^), 

S^(kj),  which  can  then  be  determined  by  setting  s = kj  in  the  integral 
equations.  A precisely  similar  situation  exists  in  "near-resonance" 
problems,  such  as  the  scattering  of  acoustic  waves  by  a long  tube,  open 
at  both  ends,  at  frequencies  near  resonance.  There  it  is  argued  that, 
although  the  function  K (t)  in  the  integral  equation  (7.33)  has  branch 
point  singularities,  the  dominant  contribution  near  resonance  comes 

from  a pole  term.  Well  away  from  resonance  it  is  anticipated  that  the  5 

dominant  contribution  to  the  Integral  comes  from  a branch  point 

singularity  which  represents  the  rather  weak  acoustic  interaction  between 

the  ends  of  the  tube.  The  functions  are  expanded  about  the  branch  point, 

and  the  Integral  term  can  then  again  be  evaluated  (approximately)  as  the 

product  of  D^(-k)  say  (where  -k  is  the  branch  point)  and  an  Integral 

which  is  expressible  in  terms  of  Whittaker  functions  (which  can  be 

further  approximated  in  most  cases) . Again  the  constant  D^(-k)  can  be 

found  by  setting  s = -k  in  the  integral  equation. 


8.  TWO-DIMENSIONAL  HALF-PLANE  PROBLEMS 


We  go  on  now  to  look  at  a simple  problem  involving  a half-plane 
embedded  in  an  acoustic  fluid.  In  the  first  place  the  fluid  will  have 
no  bulk  motion,  whereas  later  we  shall  allow  the  fluid  to  flow  over  the 
half-plane  with  uniform  subsonic  velocity,  leaving  a wake  behind  the  plate 
if  the  edge  is  a trailing  edge.  The  issue  we  want  to  examine  is  the 
following  one.  Suppose  that  the  plate  were  infinite  in  both  the  positive 
and  the  negative  x-directions  and  a wave  were  forced  to  propagate  along 
the  plate  with  some  prescribed  frequency  w and  wavenumber  q,  the  wave- 
number  q being  real  and  greater  than  the  acoustic  wavenumber  k^  = 

Thus  the  velocity  in  the  positive  y-direction  is  prescribed  in  the  form 

v(x,t)  = v^  exp(iqx  -iwt)  . (8.1) 

Then  the  solution  for  the  potential  ({)(x,y)  in  the  fluid  in  y > 0 is 
(dropping  the  factor  exp  - iwt) 

V 

4>(x,y)=  - ~ exp(iqx  - y^y)  (8.2) 

where 

Yq  = (q'  - k^)^  . (8.3) 

for  this  makes 


makes  ())-*-0asy-^+<» 

and  makes  ^ (x,0)  = v^  exp(iqx)  . 

The  field  described  by  (8.2)  is  that  of  a subsonic  trapped  surface  wave. 

No  radiation  takes  place  across  any  plane  y “ const,  because  the  pressure, 
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p “ pl(iO({i,  and  the  velocity  34i/3y  are  90  degrees  out  of  phase.  The 
energy  is  locked  in  a thin  layer  of  thickness  0(Y”')  adjacent  to  the 
surface,  and  none  escapes  as  sound. 

Suppose  now  that  the  surface  is  semi-infinite,  occupying 
(_oo<x<0,  y=0),  all  dependence  on  the  z-coordinate  parallel  to 
the  surface  edge  being  excluded.  Let  the  surface  - “>  < x < 0 still  be 
forced  to  move  with  the  prescribed  velocity 

v(x)  = v^  exp(iqx). 

Then  (8.2)  cannot  be  the  solution,  because  it  is  easy  to  see  that  since 

3(|)/3y  is  prescribed  on  y = 0,  for  some  range  of  x at  least,  (p  must  be  an 

odd  function  of  y and  since  <)>  must  be  continuous  across  the  extension 

(y“0,  0 < X < “>)  of  .he  surface,  <p  must  be  zero  there — whereas  (8.2) 

is  not  zero.  Clearly,  no  single  mode  like  (8.2),  nor  even  any  discrete 

set  of  modes  of  this  kind,  is  capable  of  making  3({)/3y  have  the  value 

v^  expiqx  on  y = 0,  x < 0 and  of  making  ^ = 0 on  y = 0,  x > 0.  The 

solution  for  ({)  must  therefore  contain  a continuous  spectrum  of  modes 

like  (8.2)  with  all  values  of  the  wavenumber.  In  particular,  it  must 

contain  modes  with  wavenumbers  a,  say,  with  a < k , and  for  these 

1 ° 

modes  the  exponential  decay  exp  {-  (q^  - k^)^y}  must  be  replaced  by 

1 ° ± 

oscillatory  behavior  exp  {i(k^  - a^)^y},  the  choice  of  + i(k^  - a*)^y 

J.  ° ° 

rather  than  - l(k^  “ oi^)^y  being  dictated  by  the  radiation  condition,  that 

the  phase  factor  exp  {i(k^  - a^)^y  - iujt)}  be  that  of  an  outgoing  wave 

° 1 

as  y -►  +“>(when  y < 0 we  take  - i(k^  - a^)^y).  Energy  is  radiated 

o 

across  a plane  y “ const,  by  such  a mode — and  we  say  that  the  energy 
which  was  trapped  in  the  subsonic  mode  (8.2)  on  an  infinite  plate  has 


51 


been  scattered  Into  other  subsonic  modes  and  also  Into  supersonic 
radiating  modes  by  the  discontinuity  In  the  surface  y » 0. 


To  analyze  the  process  of  scattering,  or  wavenumber  conversion,  by 
the  plate  edge,  we  define  a scattered  field  (p  by 

^total  ‘ 'f' 

q 

for  y 2.  0.  Because  the  derivative  the  same  value  on 

y “ 0+»  “^total  function  of  y,  and  so  It  is  enough  to 

consider  only  y 2 0*  The  scattered  field  0 is  a solution  of 

with 

-|^=0  ony=0,  x<0  (8.6) 

dy 

and  with 

V 

0 ^ exp  iqx  =0  ony=0,x>0.  (8.7) 

This  is  a typical  two-part  mixed  boundary  value  problem  which  we  may 
expect  to  solve  by  the  W-H  technique.  Two  further  conditions  are 
needed,  however,  to  get  a unique  solution  for  0.  One  comes  from  con- 
ditions expected  to  hold  as  |}c|  <»,  and  defines  the  domains  of 

analytlclty  of  half-range  transfonr:s  and  the  strip  D of  overlap.  The 
other  comes  from  conditions  at  the  plate  edge,  x ■ y » 0,  and  determines 
the  behavior  at  infinity  In  the  transform  s-plane,  and  hence  determines 
the  entire  function  arising  In  the  W-H  method.  We  shall  leave  the 
matter  of  edge  conditions  until  we  need  to  look  at  It  In  detail.  For 
the  moment  we  just  assume  that  all  functions  with  which  we  deal  have 

at  most  Integrable  singularities  at  x - y ■ 0. 
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As  to  conditions  as  x ± <»,  we  give  and  q small  positive 
imaginary  parts, 

- kj  + Ikj,  q - qj  + iq^  (8.8) 

and  then  (8.7)  gives 

(p  - 0(exp  - q^x) 

as  X -*■+",  y ■ 0.  Away  from  y * 0 we  can  expect  that  as  x -►+®  p will  take 
the  form  of  an  outgoing  cylindrical  wave, 

1 

((i  - r ^ exp(lk^r)  f(6)  - 0(exp  - k^x)  (8.9) 

It  then  follows  that  all@  functions  linearly  related  to  ()>  and  its 
y-derivatives  will  be  analytic  in 

: Im  s > - min(q2.kj)  (8.10) 

As  X -+--00,  contains  an  exponentially  growing  part  which  we  have 

split  off  in  (8.4),  so  that  4)  should  behave  like  an  outgoing  cylindrical 
wave, 

P - 0(exp  - kjxl)  (8.11) 

as  X Then  all Q functions  will  be  analytic  (and  with  algebraic 

behavior  at  infinity)  in 

R_  : Im  s < + kj,  (8.12) 

and  the  strip  D is 

D:  - min(qj,,k2)  < Im  s < + k^ . (8.13) 
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For  y > 0 we  can  apply  Fourier  transforms  to  (8.5),  and  can 
Integrate  the  3^4i/dx^  exp  Isx  term  by  parts  twice  with  no  contribution 


from  X * ± ^ provided  s lies  In  D.  This  gives 


where 


Y » (s^  - , 

s o 


(8.14) 


(8.15) 


Introducing  the  square  root  function  whose  behavior  In  the  complex 
s-plane  holds  the  key  to  many  aspects  of  acoustic  diffraction  and 
scattering  processes. 

± 1 

The  function  Y = (s  - k )^(s  +k  )^  has  branch  points  at  s = ± k , 
s o o o 

and  branch  cuts  must  emanate  from  these  points  to  form  a barrier  which 

must  not  be  crossed.  We  can  either  make  a cut  from  + k to  -k  , or 

o o 

we  can  make  a cut  from  + k to  (in  any  direction)  and  a cut  from 

° 1 1 

-k  to  (in  any  direction).  If  the  values  of  (s  + k )^,  (s  - k )^ 

O 0 0 

are  specified  at  any  point  in  the  s-plane  (not  necessarily  the  same 

point  for  the  two  functions)  and  the  branch  cuts  are  fixed,  then  a 

1 _i_ 

unique  value  of  (s  + k )^,  (s  - k )^  is  obtained  by  starting  at  the  given 

o o 

point  and  moving  to  any  desired  point  without  crossing  any  branch  cut, 
and  Insisting  that  the  function  change  continuously  from  its  initial 
value.  Figure  4 gives  various  possible  choices  of  branch  cuts.  In 


addition  to  the  choice  of  branch  cuts  we  shall  take  each  of  (s  ± k )* 

I 1 ° 

to  be  the  branch  which  behaves  like  + s^  (rather  than  - s^)  when  s is 
large  and  positive. 

Now  in  our  problem  we  know  that  <l>  must  be  analytic  in  D,  and  since 


the  general  solution  of  (8.14)  is 


1 


♦(s.y)  “ A(s)  e ^8^  + B(8) 


(8.16) 


that  will  only  be  possible  If  the  branch  cuts  from  ± do  not  enter 

the  strip  D.  Thus  the  cut  from  +k^  must  go  to  Infinity  above  the  strip, 

that  from  -k^  to  Infinity  below  the  strip.  No  further  specification  of 

the  cuts  need  be  made  at  this  stage,  because  we  shall  find  a solution 

for  s In  D,  and  the  values  of  y are  already  fixed  for  s In  D by  the 

1 J.  ® 

requirement  (s  ± k^)^  ~ + s^  as  s -♦•  + “>  and  by  the  general  location 
of  the  cuts. 

We  can  now  see  that 


0 < arg(s  + k^)  < tt 


-TT  < arg(s  - k^)  < 0 


for  all  s In  D,  and  therefore 


(8.17) 


- -j  < arg(s^ 


< I 

o 


or  equivalently 


Re  Y >0  for  all  s In  D 

E 

which  Is  the  essential  property  H Is  possible  to  choose  the 


(8.18) 


branch  cuts  so  that  (8.18)  holds  throughout  the  entire  complex  (cut) 
plane,  but  there  Is  no  need  for  this  since  at  the  moment  we  are  concerned 
only  with  values  of  s In  D.  Then  the  only  possible  form  for  B(8)  In 
(8.16)  Is  B(s)  ■ 0,  otherwise  4>  would  be  Infinite  as  y -*■  + “.  Thus 


<K8,y)  - A(s)  exp(-Y  y) 


(8.19) 
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and  Che  original  radiation  condition  of  outgoing  waves  is  seen  to  be 
equivalent,  for  > 0 and  s in  D,  to  the  condition  that  4(s,y)  •+•  0 as 
y -*■  + *. 


The  boundary  condition  (8.6)  involves  4'(s,0),  where  the  ' indicates 
3/dy,  and  therefore  we  differentiate  (8.19)  to  get  the  pair  of  equations 
(equivalent  to  the  differential  equation  plus  the  radiation  condition) 


4>(s,0)  * <I>^(s,0)  + <t_(s,0)  « A(8) 

♦ '(s,0)  - 4>'(s.O)  + *'(8,0)  - - Y A(8) 

•T  - S 


(8.20) 


where,  for  example. 


.0)  e dx  , 


♦ (8,0)  - ({i(x,( 

“ OD 

♦;(s,0)  - J"  1^  (X,0)  e^®*  dx. 

In  (8.20)  two  of  the  functions  aie  known.  From  (8.6)  we  have 


etc. 


<*»'(s,0)  ■ 0 (for  8 in  R ) 
while  from  (8.7)  we  have 

v«.'»  ■ 

for  Im  s > - qj,  l.e.,  for  seR^^.  Eliminating  A(s)  between  the  two 
equations  in  (8.20)  and  using  (8.21)  and  (8.22)  gives 

- iv 


(8.21) 


(8.22) 


K(s)  ♦;(s,o)  + ♦_(s,o)  - - -3% 


(8.23) 
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a standard  form  of  W-H  equation,  with 


K(8)  - - (s*  - 

'a 


(8.24) 


The  factorization  of  K Into  factors  analytic  and  non-zero  in  R 
and  of  algebraic  behavior  at  Infinity  there  is  again  obvious: 


K (s)  - (s  + k ) ^ , K (s)  - (s-  k ) , 

+ o - o 


and  after  division  of  (8.23)  by  K_(s)  we  again  have  the  additive  split 
of  the  function  l/(s  + q)  K (s)  to  make.  We  thus  arrive  at  the 


equation 


Vs)  Yq  K_(-q)(8~"qy 

<^-(8,0)  iv^  j ^ 

K (s)  Y (s  + q Ik  (s)  K (-q) 

q \ - 


(8.25) 


* an  entire  function  E(s). 

We  anticipate  that  E(s)  will  be  a polynomial,  and  consider  the  implica- 
tions that  the  degree  of  the  polynomial  be  N,  E(s)  ■■  a^s^  + a^s*^  +...+ajj. 

iv  (s  + k )^  i.  « 

((,•  (8,0)  T + (s  + k )^  {a  s"+. . .} 

y K (-q)  (s  + q)  o o 

q - 

.1  i. 

- 0(8  *)  + 0(a  S*”””) 

O 


as  s “ in  R^.  The  first  term  vanishes  algebraically  at  “ and  is 
therefore  the  half-range  transform  of  a function  with  at  most  an  Integrable 
singularity  at  x - (H-.  In  fact  it  follows  from  Watson's  lemma  that  the 
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_1_ 

corresponding  value  of  3(|)/3y  on  y • 0,  x > 0 is  0(x  as  x -►  0+. 

The  second  term  grows  as  s and  must  arise  as  the  generalized  Fourier 

transform  of  a function  which  has  a non-integrable  singularity  at  x * 0. 
According  to  Llghthlll  [ll,p.43] 


r 


X ^ exp  Isx  dx  = e 


T I iX  1 

(-X)! 


for  real  s,  and  the  correct  Interpretation  of  this  in  is 

T A-, 

e (-X)!  s ‘ 

where  the  branch  cut  for  the  function  s^~*  is  to  go  from  0 to 
in  the  lower  half-plane. 


(8.26) 


Hence 


_i^ 

♦^Cs.o)  - 


<=>  (x,0)  = 0(x  ^ as  x -►  0+. 

3y 

—>2 

Therefore  the  velocity  has  a singularity  at  least  as  bad  as  x near 
x = 0,  and  the  kinetic  energy  in  a small  region  around  x ■ 0 will 
diverge  to  infinity. 

We  argue  that  this  singularity  is  unacceptable,  and  choose  the 
solution  corresponding  to 

E(s)  S 0 , 

\ 

thtts  giving  the  least  singular  behavior — like  x — in  the  velocity  at 
X * 0.  Note  that  it  is  impossible  to  Impose  a Kutta  condition,  that 
the  velocity  be  finite  (except  by  abandoning  the  radiation  condition, 
or  allowing  to  be  discontinuous,  and  that  cannot  be  permitted  in 
static  fluid) . To  see  what  kind  of  a pressure  field  exists  near  x ■ 0 
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we  have 


iv 


iv 


* (8,0) 


Yq(s  + q) 


,1  - 


(s  -k^)"  K_(-q) 


On  y =*  0,  X > 0 we  know  from  (8.7)  what  <|)  should  be,  and  this  can  be 
confirmed  from  the  expression  for  On  y “ 0,  x < 0 we  close  the  Inverse 

Fourier  integral  path  in  the  upper  half-plane.  The  pole  s = -q  lies 
outside  the  contour  and  makes  no  contribution,  so  chat  we  only  need  to 
examine  the  second  contribution  to  <I>_  as  s in  R . For  this  contri- 

bution 

_ 3 

it  (s,0)  ~ s ^ 


and  so 


(()(x,0)  ~ (-x) 


(times  some  coefficient)  as  x -►  0-.  Thus  the  pressure  and  the  pressure 

1 

jump  both  vanish  like  (-x)  near  the  plate  edge.  Note  that  although 
the  pressure  jump  does  vanish  at  the  plate  edge  (which  would  be  regarded 
in  aerofoil  theory  as  the  satisfaction  of  a Kutta  condition) , the 
velocity  is  nonetheless  infinite  at  the  edge. 

In  summary,  the  least  singular  solution  has 
1 \ 


4.  - 0(x'^) 

1 

V4)  - 0(x ’’') 


(8.27) 


near  the  edge,  and  conditions  of  this  kind  are  often  Imposed  at  the 
outset  as  edge  conditions.  It  seems  preferable  not  to  anticipate  the 
edge  behavior  in  advance,  but  to  follow  the  W-H  method  through  as  far 
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as  (8.25),  and  Chen  to  see  In  each  particular  case  what  behavior  muse 
hold  near  Che  edge  and  vhaC  freedom  exists  for  minimizing  singular 
behavior  (as  we  shall  do  In  ^9).  In  some  cases.  In  particular  In  recent 
work  by  Rawlins  [IS]  on  diffraction  of  an  acoustic  wave  by  a half-plane 
which  Is  "sound-hard"  on  one  side  and  "sound-absorbing"  on  the  other, 
the  edge  conditions  are  not  at  all  obvious,  being  In  fact 

1 

4>  = O(x^) 

_ 3 

74>  « 0(x 

To  justify  acceptance  of  a solution  with  certain  edge  conditions 
one  has  to  go  beyond  the  simple  linear  Invlscld  wave  equation  used  here, 
or  beyond  Che  simple  zero  thickness  model  of  Che  boundary.  For  example, 
one  can  look  at  linear  acoustic  propagation  with  viscous  effects 
Included,  or  one  can  look  at  Invlscld  propagation  around  a surface  which 
Is  Chin  compared  with  any  ocher  relevant  length  scale  but  which  has  a 
smoothly  rounded  edge.  IL  can  be  proven  (though  the  proof  has  not  yet 
been  published)  that  our  solution  with  conditions  (8.27)  Is  the  unique 
one  which  can  be  matched  to  an  "Inner  solution"  In  which  either  viscous 
forces  or  the  continuous  curvature  of  the  boundary  lead  to  finite 
velocities  everywhere.  That  Is  a rather  special  kind  of  proof,  however, 
and  we  shall  refer  In  f>9  to  Che  unavailability  of  a comparable  proof 
when  there  Is  uniform  subsonic  flow  past  the  radiating  half-plane. 

This  discussion  of  edge-conditions  completes  the  formal  determina- 


tion of  Che  field  as 


where  c runs  from  - “ to  + «>  in  D,  l.e.,  above  s » -q,  above  the  branch 


point  8 “ -k  and  below  the  branch  point  s = +k  . This  holds  for  y > 0, 
*0  o 

while  for  y < 0 we  use  “ “‘^total^^^’  integral  here  can 

actually  be  evaluated  in  closed  form,  in  terms  of  Fresnel  integrals, 
though  the  details  are  complicated,  and  are  nothing  to  do  with  the 
W-H  method.  We  refer  the  reader  to  [5,10]  for  descriptions  of  the 


details,  and  remark  simply  that  the  distant  radiating  acoustic  field 
can  be  estimated  asymptotically  from  the  formula 


F(s)  exp(-isx  - Y^y)  ds 


(8.29) 


1 

/2k  71 

_ i j exp(ik^r  - 7Ti/4)  sin6  F(~k^  cosS) 


in  which  x - r cos0,  y = r sin0,  0 <.  0 i tt.  (This  formula  may  give 

apparent  infinities  in  particular  angular  directions,  and  near  those 

angles  it  is  necessary  Lo  use  more  refined  approximations.)  To  apply 

1 

it  here  we  have  to  Interpret  (-k  cos6  -k  )^,  which  we  take  as 
, o o 

-ik  ^(1  + cos0)^  because  the  arg  of  -k  cos0  - k is  equal  to  -tt  for 
o o o 

all  0 between  0 and  1T.  We  also  need 

■ <-■'  - • 
again  because  arg  (-q  - k^)  • -tt  . Then  we  have 


1 

,2 


1 ik  r -771/4  sln0/2 


(q  - k^)T 


(q  - k^  COS0) 


which  gives  the  level  and  directivity  of  the  scattered  acoustic  field. 

In  the  next  section  we  look  at  the  same  problem,  but  with  sub- 
sonic mean  flow  over  the  surface. 


(8.70) 


9.  HALF-PLANE  PROBLEMS  WITH  MEAN  FLOW;  WAKES  AND  THE  KUTTA  CONDITION 

Consider  again  the  half-plane  y » 0,  x < 0 with  prescribed  velocity 
v(x)  ■ exp(lqx),  but  suppose  now  that  there  is  uniform  parallel  flow 
at  the  same  speed  U on  both  sides  of  the  plate,  the  plate  edge  being  a 
trailing  edge.  The  potential  ^he  convected  wave  equation 

- ‘I  l7‘)  ♦total  ■ “ 

and  the  boundary  condition  on  the  plate  is  that  of  continuity  of 
displacement  (not  of  normal  velocity) , giving 


(-♦“  * “ ii) 


n “ T" 

3y 


total 


where  n “ v/(-lu))  is  the  surface  displacement, 


Thus 


at)! 

-lo)  -5—' 
3y 


total 


(-iu)  + iuq)  v^e 


iqx 


on  y “ 0,  X < 0. 


(9.2) 


Because  of  the  presence  of  mean  flow  there  is  a new  possibility  on  the 

extension  of  the  plate.  A wake  can  exist  there,  across  which  there  may 

be  a discontinuity  in  tangential  velocity  as  long  as  the 

normal  displacement  (and  hence  here  also  the  normal  velocity)  and  the 

pressure  are  continuous  across  the  wake.  Take  a single  Fourier  com- 
iXx 

ponent  Ae  of  ^ x > 0;  then  for  y ■ 0-,  x > 0, 

iXx 

below  the  wake,  the  corresponding  component  is  -Ae  , because 

be  an  odd  function  of  y.  The  pressure  jump  across  the  wake  associated 

with  this  particular  component  of  potential  is 


p(x,(H-)  - p(x,0-)-(-p)(-iu)  + U (2Ae^^’‘) 
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and  this  vanishes  either  if  A - 0 (vrtiich  implies  continuity  of 

and  hence  absence  of  the  wake)  or,  for  any  value  of  A,  if  X ■ u)/U.  Thus 

the  general  condition  on  across  the  wake  is  that 

iU)x/U 


'Ototal 


(9.3) 


for  some  value  of  A.  The  tangential  velocity  is  (lwA/u)e^*^'^  above 

the  wake,  (-la)A/U)e^*^^^'^  below  it,  and  hence  the  wake  is  an  oscillatory 

vortex  sheet  of  strength  2iuA/U,  modulated  by  the  phase  factor 

exp  ^i  ^(x  - Ut)J  , which  shows  that  any  element  of  the  vortex  sheet 

propagates  downstream  at  the  flow  speed  U. 

We  wish  to  solve  the  two-part  mixed  boundary  value  problem  posed 

by  (9.1-9. 3),  subject  here  to  the  restriction  to  subsonic  flow 

(M  - U/c  < 1)  and  to  a radiation  condition,  but  leaving  open  the  issue 
o 

of  edge  conditions.  We  first  write 


•'total  ^(-q) 
where  we  shall  write  generally 


D(-q)  exp(iqx  -5(-q)y)  + (j) 


(9.4) 


s*  - (k  + Ms)2 
o 


1 

T 


(9.5) 


^“9)  - Cs  . .q 

D . 1+ 1+^ 

s 0)  k 

o 


D(-q)  - 1 


. SH 


U) 


(9.6) 


The  function  ^ replaces  the  no-flow  case,  while  is  a kind 

of  Doppler  factor  for  wavenumber  s.  The  object  of  writing  in 


the  form  (9.4)  is  that  the  incident  field  associated  with  the  velocity 


on  an  infinite  plate  has  been  split  off;  that  field  is  exponen- 
tially large  as  x and  the  ({)  that  remains  is  assumed  to  be  just  an 

outward  propagating  wave. 

We  now  have  the  following  problem  for  <}>! 


(-1.  + U 


(9.7) 


= 0 (y  = 0,  X < 0)  (9.8) 

dy 

V D(-q) 

0 = A exp(^^)  + — exp(iqx)  (y  ® 0,  x > 0)  (9.9) 


and  we  have  to  examine  conditions  as  x ± Write  k = k,  + ik, , 

o 1 2 

u)/U  = (kj  + ik^)/M,  q = q^  + As  x (Ji  must  behave  like  an  out- 

going acoustic  wave,  whose  decay  will  be  like  that  of  a plane  wave 
propagating  against  the  flow,  so  that 


(9.10) 


As  X ->■  + «“,  ())  may  behave  either  like  an  acoustic  wave  propagating  with 
the  flow 


(J)  ~ exp 


iu)t  + 


ioix  I 

c + U j 
o ' 


~ exp 


k^x 
1 + M 


(9.11) 


or  may  behave  as  it  does  in  the  wake 


(9.12) 


fk^x 

^ - «P  - j — 

or  In  the  way  associated  with  the  forcing  field  exp  Iqx, 

<|>  ~ exp  {-  x} 

It  follows  from  these  chat 


R Is  Im  s < + — 

"•  1 “ M 

and  that  the  strip  Is 

I 

provided  M < 1. 

Take  Fourier  transforms  of  (9.7)  for  y > 0 to  get 

<I>(s,y)  = A(s)  exp  (-^  y)  + B(s)  exp(+^  y) 

s s 


Im  s < + 


1 - M 


/ k, 

R Is  Im  s > - min  I - ■ 
+ \ 1 + M 


(9.13) 


(9.14) 


(9.15) 


C Is  defined  In  (9.5).  We  take  a cut  from  s ■ + k /(I  - M)  to  Infinity 
s o 

above  the  strip  D and  one  from  s » -k  /(I  + H)  to  Infinity  below  D,  and 

± ° ±1 
define  [s  - (k  +M8)]^  to  be  the  branch  which  behaves  like  + (1  - M)^  s^ 

° 1 
when  s Is  large,  real,  and  positive,  [s  + (k  +M8)]^  to  be  the  branch 

1 1 o 

behaving  like  + (1  + M)^  s^  for  large  real  positive  s.  With  these  branch 

cuts  and  choices  of  branch.  Re  r >0  for  all  s In  0,  lust  as  for  y 

s s 

In  ^8.  Then  since  (9.15)  only  holds  In  D we  have  to  have  B(8)  ■ 0, 
and  now  we  can  differentiate  (9.15)  with  respect  to  y,  put  y • 0-t- 
and  eliminate  the  function  A(8) . Thus 
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♦{.(s,(H)  + 'l>^(8,0+) 


- Cgt4>^(8,(H)  + 4>_(8,Of)] 


From  (9.8)  we  have  <t''(8,0+)  ■ 0 while  from  (9.9)  we  have 


lA 

♦+(8.0+)  - ^ ^ + C(-q)(8  + q) 


and  hence  we  get  a W-H  equation 


(t;(8,0+)  + $_(s.(H-)  + +^k^/M)  5(-q)(s  + q) 

We  define  the  factore  4^(s)  C (s)  by 


^^(s)  ••  [s  + (k^  + Ms)]^ 


C (8)  - [8  - (k  + M8)j5' 
- o 


(9.16) 


(9.17) 


and  by  the  usual  route  arrive  at 

k ^'^o 

♦^(8,0+)  + (8  + ,,^/M)  iT^  f;(-q)(s  + q) 


-i(8)  l(s.o^)  - Csfr/M) 


Iv  D(-q) 

s(4)(.~7) 


(9.18) 


■ an  entire  function  E(s),  which  must  be  a polynomial  because  of  the 
algebraic  behavior  of  all  functions  Involved  here.  In  fact  we  must  have 
E(8)  - 0,  for  If  E(8)  were  even  a constant,  the(5p8>^t  of  (9.18)  would 
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4^(8, Of)  * 0(8*)  at  infinity 

which  corresponds  to  the  strong  singularity  V(])  ^ 0(x  *).  With  that 

choice  of  E(s)  the  solution  is  still  not  unique,  because  the  wake 

strength  A is  undetermined.  We  can  argue  that  conditions  are  such 

that  there  should  be  no  jump  in  tangential  velocity  across  the  extension 

i 

of  the  plate,  in  which  case  A • 0,  V(fi  - 0(x  ^)  and  (|)  = 0(x*)  near  the 
plate  edge  and  we  have  a situation  essentially  the  same  as  was  examined 
in  ^8.  Because  the  radiation  is  emitted  into  uniformly  moving  fluid 
the  "stationary  phase"  formula  (8.29)  is  not  immediately  applicable,  but 
we  shall  show  in  a moment  how  it  can  be  generalized  to  the  moving  fluid 
case. 


We  can  alternatively  argue  that  a wake  will  adjust  Itself  to 
eliminate  the  high  velocities  which  would  otherwise  exist  at  the  trail- 
ing edge,  and  that  a Kutta  condition  of  flnlteness  of  the  velocities 
at  the  edge  should  be  imposed  whenever  possible.  The  physical  basis 
for  such  a condition  in  unsteady  trailing  edge  flow  is  a matter  of 
controversy  at  the  moment,  but  that  does  not  concern  us  here  as  our 
Interest  is  merely  in  seeing  1£  and  how  such  a condition  can  be  applied 
in  this  model  problem. 

Expand  the  ©part  of  (9.18)  as  s -*•  “ in  R^.  We  have 


$'^(s,0)  ~ ^1AC_(-  -^) 


iv^D(-q) 


C(-q) 


(-q)  j 


(1  + M)*  s 


± i 
2 »-  T 


+ 0(8"'^) 


(9 


the  term  given  explicitly  corresponding  to  V4)  * 0(x  ^)  the  second  to 
1 

V(|>  * O(x^).  Thus  we  can  impose  a Kutta  condition,  that  the  velocities 
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.19) 


be  finite  at  the  edge  (note  that  the  velocities  in  the  part  of 

split  off  in  (9.4)  are  finite,  but  non-zero,  near  the  edge)  by  choosing 

a wake  strength 


D(-q) 

l(-  t)  «+‘-‘ 


(9.20) 


It  is  easy  to  see  that  the Q part  of  (9.18)  contains  terms  like  (s  + q)”*, 

(s  + k /M)”*  which  (as  in  68)  make  no  contribution  to  ({)  for  x < 0, 
o 

y ■■  O-t*  and  that  with  the  choice  (9.20)  the  first  term  in  the  expansion 

_ 5 

of  i (s,0+)  as  8 " in  R is  0(s  ^) . This  corresponds  to 


4)  “ O(x^) 

“ /a  3 \ 

the  edge,  though  the  pressure  is  O(x^)  because  p “ ^ 

The  expression  for  the  potential  in  y > 0 is  found  to  be 


j F(s)  exp(-isx  - C^y)  ds 
^ c 


k iv  D(-q) 
f / _Ov  . o 

C+(-q)(s+q) 


(9.21) 


In  the  exponential  factor  we  write 


-isx  - [s^  - (k  + Ms)^]^  y 


- isx 


SO  that  if  we  define 
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1 

4>^(8,0+)  ■ 0(8^)  at  Inflntty 

j_ 

which  corre8pond8  to  the  8trong  8lngularlty  V()>  > 0(x  ^).  With  that 
choice  of  E(8)  the  aolutlon  la  atlll  not  unique,  because  the  wake 
strength  A Is  undetermined.  We  can  argue  that  conditions  are  such 

that  there  should  be  no  Jump  In  tangential  velocity  across  the  extension 

1 1 
— "S’ 

of  the  plate,  In  which  case  A"0,  V(j)»0(x  ) and  ()>  “ 0(x  ) near  the 
plate  edge  and  we  have  a situation  essentially  the  same  as  was  examined 
In  ^8.  Because  the  radiation  Is  emitted  Into  uniformly  moving  fluid 
the  "stationary  phase"  formula  (8.29)  Is  not  immediately  applicable,  but 
we  shall  show  In  a moment  how  It  can  be  generalized  to  the  moving  fluid 
case. 


We  can  alternatively  argue  that  a wake  will  adjust  Itself  to 
eliminate  the  high  velocities  which  would  otherwise  exist  at  the  trail- 
ing edge,  and  that  a Kutta  condition  of  flnlteness  of  the  velocities 
at  the  edge  should  be  imposed  whenever  possible.  The  physical  basis 
for  such  a condition  In  unsteady  trailing  edge  flow  Is  a matter  of 
controversy  at  the  moment,  but  that  does  not  concern  us  here  as  our 
Interest  Is  merely  In  seeing  1£  and  how  such  a condition  can  be  applied 
In  this  model  problem. 

Expand  the®  part  of  (9.18)  as  s In  R^.  We  have 


4>' 


^(8,0)  - ^1A5_(-  ^)  + 


Iv  D(-q)  C (-1 


e(-q) 


(-q)^ 


i.  I 
(1  + M)*  8* 


+ 0(8“'^) 


(9 


1 

-Tx 


the  term  given  explicitly  corresponding  to  V^)  > 0(x  ^)  the  second  to 
1 

74)  ■ O(x^).  Thus  we  can  Impose  a Kutta  condition,  that  the  velocities 


.19) 


67 


a * s 


K * T1J2 
o 1-M 


* 

X = R cos@,  (1  - y = R sin@ 

then 

exp  |-io  R COS0-  (o^  - R sin0jda 

which  is  now  capable  of  asymptotic  estimation  by  the  formula  (8.29)  for 
static  fluid. 

To  make  the  algebra  less  complicated,  suppose  the  Mach  number  M 
is  small  so  that  can  be  neglected  compared  with  unity.  Then  the  field 
associated  with  the  second  contribution  in  (9.21)  turns  out  to  be 


iv^(l  - Mq/k^) 


(<,  - 'lo_J 

\ l+M/ 


sin0/2  exp  (ik  r - ik  Mx  - (9.23) 


fq  - k (cos0  - M)  ] 
o 


there  being  no  difference  between  r and  R or  between  0 and  0 if 
« 1.  The  field  (9.23)  is  a trivial  modification  of  (8.30). 

The  distant  field  associated  with  the  first  term  in  (9.21),  with 
the  Kutta  condition  value  (9.20)  for  A comes  out  as 


-(y 


iv  (1  - Mq/k  ) 


(1-  M cose) 


8ln0/2  exp{ik  r - ik^  Mx  - ^ } 
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and  the  ratio  of  the  wake  generated  field  (9.24)  to  the  field  (9.23) 
which  would  exist  In  the  absence  of  the  wake  Is  (since  q > k ) roughly 


(9.25) 


This  shows  that  If  the  plate  wave  travels  at  speed  co/q  between  the  flow 

speed  and  the  sound  speed,  the  wake  contribution  to  the  far-fleld  Is 

negligible.  If,  on  the  other  hand,  the  plate  wave  travels  more  slowly 

than  the  flow  speed,  Mq/k^  exceeds  unity  and  the  sound  field  when  the 

Kutta  condition  Is  Imposed  exceeds  that  In  the  absence  of  any  wake. 

This  Is  not  to  be  regarded  as  a general  conclusion,  for  In  other 

problems  the  wake  sound  very  nearly  cancels  the  primary  edge  field.  That 

In  fact  happens  here,  for  when  q becomes  close  to  k /M  It  can  be  seen 

o 

that  the  fields  (9.23)  and  (9.24)  have  small  but  equal  and  opposite 
values.  Whether  the  extra  wake  field  dominates,  or  mainly  cancels 
the  primary  edge  sound  field  thus  depends  very  much  on  details  of  the 
basic  excitation. 
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10.  CONSTRUCTION  OF  W-H  SPLIT  FUNCTIONS 


In  this  section  we  first  outline  a general  method  for  effecting 
either  the  additive  or  the  multiplicative  decomposition  of  a function 
analytic  In  a strip,  and  then  ve  set  down  some  properties  of  the  functions 
which  arise  frequently  In  acoustics  problems.  Finally,  we  record  the 
corresponding  properties  for  strictly  Incompressible  flow  problems. 

A.  Cauchy  Integrals 

Let  F(s)  be  analytic  in  some  strip  D;  e < Im  s < e and  let  R 

1 2 + 

denote  the  domain  Im  s > e^,  R_  the  domain  Im  s < e^.  Suppose  also  that 
|F(s)1  -*•  0 uniformly  as  fast  as  |s|~^for  some  X > 0 as  |s|  -*■<*>  within 
any  closed  region  within  D,  l.e.,  as  |s|  ->■<»  with 


Ej  < Oj  <.  Im  s £ Oj 


< e. 


Then 


F(a)  - F^(s)  + F_(8) 


for  s In  D,  where 


Is  analytic  and  bounded  In  R^, 


F (s) 


t litl 


2711 


I 


t -s 


dt 


(10.1) 


(10.2) 


(10.3) 


Is  analytic  and  bounded  in  R . 

The  path  runs  from  - * to  + » In  D below  t ■ s,  tihile 

— /'^^runs  from  - <*>  to  + * In  D above  t ■ s. 

Suppose  further  now  that  |F(s)|  -*-0  uniformly  as  fast  as  |s| 


for  some  p > 0 as  |s|  -*-<>i>ln  the  strip  D.  Then  | F(t)dt  converges 
absolutely,  and 


I 
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nt)dt-0(,-) 


(10.4) 


as  |8|  -►“•In  R_j_,  while 


F (s)  . 


-A.f 

2Tris  J 


F(t)  dt  = 0(s~‘) 


as  I s I -►  “>  in  R . 

Without  going  Into  fine  details,  (10.1)  is  proved  by  applying 
Cauchy's  theorem 

to  a contour  lying  within  D and  enclosing  the  point  s in  D (F(s)  being 
defined  only  for  s in  D in  the  first  Instance).  The  contour  is  then 

deformed  to  consist  of  a rectangle  with  — ► as  its  lower  side 

•4 — ^ as  its  upper,  and  with  the  ends  of  these  sides  joined  at 

infinity  by  short  sides  parallel  to  the  Imaginary  axis.  In  the  limit 
these  short  sides  (of  finite  length,  less  than  £2  ~ make  no 

contribution  to  the  Integral,  so  that 


v(s)  . f litl  JL  f F(t)_ 

2Tri  J t - s 27ri  J t - s 


where  now  both  paths  of  integration  run  from  - » to  + «>.  But  then, 
according  to  the  basic  theorem  of  complex  variable  analysis,  the  first 
term  defines  an  analytic  function  as  s varies  without  crossing  the 
integration  path,  l.e.,  it  defines  an  analytic  function  in  the  upper 
half-plane  above  the  Integration  path.  Similarly,  the  second  term 
defines  a function  analytic  everywhere  below  the  integration  path 


For  the  behavior  at  infinity  we  have 


(10.5) 
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Without  any  real  loss  of  generality  we  can  take  the  path  of 
to  be  the  real  axis,  and  s^  ^ 0 in  so  that 

I f tP(t)  I < (*“  Itl  lF(t)l  <lt 

' L At  - s.)‘ 

Let  |s|  “ in  along  a ray  with  s^  = Ks^,  and  divide  the 

integration  at  points t “ ±M.  M is  chosen  so  that  |s|  » M, 
that  for  |t|  >.  M,  |F(t)|  < c|t|  ^ ^ for  some  constant  c and 
Then  on  (-M,+M) 

1“  |t|  lF(t)ldt  - Ip  |j|  |p(t)| 

-M  At  - Kspt  + s' 

and  the  integral  is  finite  and  independent  of  s,  while  on  say  we 

put  t * Ks^  + s^tanS  to  get 


Integration 


range  of 
but  so 
some  A > 0. 


r 


t|F(t)ldt  _^c_ 

/(t  - KSj,)^  + 8\  ®2 


r 


K> 


sec  6 

(K  + tan0)^ 


d0 


and  again  the  Integral  is  convergent  and  Independent  of  s.  Thus  under 
these  conditions 

F^(s)  “ 0(s“‘)  at  Infinity  in  R^. 

The  above  is  hardly  a proof,  but  it  can  be  rlgorlzed.  In  any 
application  the  behavior  at  infinity  should  be  checked  out  carefully  in 
each  case. 
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The  Cauchy  Integral  formulas  (10.2,  10.3)  enable  the  product 
decomposition  K(s)  * K^Cs)  K (s)  to  be  effected  by  taking  logarithms. 

Suppose  K(s)  is  analytic  in  the  strip,  and  that  |K(s)|  1 uniformly 

as  |s|  -*•  “ in  the  strip.  Suppose  further  that  K(s)i<0  in  the  strip. 

Then  F(s)  ■ f-n  K(s)  is  analytic  in  the  strip  for  any  branch  of  the 
logarithm,  and  can  be  decomposed  as  in  (10.1).  Define 

K^(s)  - exp  F^(8),  K (s)  ■ exp  F_(s)  . (10.6) 

Then 

K(s)  - K^(s)  K_(s)  (10.7) 

for  s in  D,  and  K^(8)  are  analytic  and  non-zero  throughout  R^, respectively. 

This  decomposition  is  unique  up  to  multiplication  of  say  K^(s) 
by  a non-zero  entire  function  and  division  of  K (s)  by  the  same  function. 

It  may  be  necessary  to  use  this  freedom  to  remove  non-algebralc  behavior 
at  infinity  of  K^(s)  in  certain  cases.  Noble  [5]  gives  several  examples 
of  this  minor  difficulty. 

When  K(s)  is  even,  the  factors  K^(s)  as  defined  by  (10.6),  (10.2), 
and  (10.3)  have  the  property 

K^(-8)  - K_(8)  (10.8) 

If  the  split  is  achieved  by  some  way  other  than  use  of  Cauchy  Integrals 

it  may  be  necessary  to  adjust  the  functions  before  (10.8)  holds.  Etor 

1 

example,  if  K(s)  ■ (s*  - k*)^  then  the  "obvious"  split  is 

1 1 
K^(8)  - (8  +k^)'^,  K_(s)  - (s  - k^)"*"  , 

but  K,(-8)  - (s  - k )^,  so  that  we  need  to  redefine  K. (s)  as 

+ o z 


7A 


1 


K (s)  - (8  + k )^  K (s) 

+ o - 

in  order  that  (10.8)  will  be  satisfied. 


Tri/4 

e 


B.  Decompositions  Related  to  the  Square  Root  Function 

As  Just  noted,  the  multiplicative  decomposition  of  y = (s*  - 

s o 

into  factors  analytic,  non-zero  and  algebraic  at  infinity  is 

1 

K. (s)  * a(8  + k 

° 1 

K (s)  = a"^(s  - k )*  (10.9) 

- o 

for  any  constant  a.  If  it  is  useful  to  require  that  K^(-8)  > K (s) 
then  a should  be  chosen  as  e 

For  the  convected  wave  equation  (with  subsonic  convection  veloci- 
ties) Yg  is  replaced  by 

= [s^  - (k  + Ms)^]^ 

o O 

for  which  the  multiplicative  split  is 

1 k 1 
K^(8)  - a(l  - M*)"^  (8  + -j^)^ 

k J. 

K_(8)  - a-‘(8  - ^)*  (10.10) 

for  any  constant  a.  Here  of  course  we  cannot  make  K^(-s)  ■ K (s). 

Because  the  factors  K^(8)  K_(s)  for  K(8)  are  analytic  and  non-zero, 
the  split  for  1/K(s)  is  given  by  (l/K^(8)) , (1/K_(8) ) . 

The  additive  decomposition  of  and  functions  related  to  it  arises 
very  frequently.  Noble  [5]  gives  several  ways  of  calculating  the  split 
functions  and  several  representations  of  those  functions.  Here  we  will 
Just  verify  that  if 
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r 


1 


then 


^+(s)  “ ~ 

p (s)  ■ — cos  ' (-s/k  ) 

- TT  O 


P+(s)  + P_(s)  - Y. 


(10.11) 


(10.12) 


for  s in  D,  and  P^(8)  are  analytic  in  with  a certain  behavior  at 

infinity  which  will  be  determined  shortly. 

Firstly  it  is  necessary  to  define  the  function  cos”^  (s/k  ) for 

o 

complex  s.  If  s/k^  is  real,  cos~*  (s/k^)  is  defined  as  the  branch  for 

which  cos”'  (s/k  ) “ v/2  when  s/k  «■  0,  i.e.,  cos”'  (s/k  ) lies  between 
o o o 

0 and  IT  when  s/k  is  real  and  between  ±1.  Let  6 * cos”'  (s/k  ).  Then 
o o 


s/k 


i3  . -iB 

e + e 


and  so 


iB  = Jin 


Jin 


s (s'  ” k')' 
^ + 2_ 


(^) 


where  Yg  ” (s^  ~ with  the  branch  cuts  as  already  discussed.  The 

logarithm  here  is  defined  to  have  its  principal  value,  i.e.,  in  z is 
such  that 

-Tf  < Im  Jin  z<  -Mr  (10.13) 

with  a branch  cut  along  the  negative  real  z-axls.  Now  take  s ■>  0;  the 

TT  i 

corresponding  value  of  Y is  "ik  and  hence  iB  * Jln{±(-i)}  * + — » so 

o 2 

that  B * tt/2  if  we  choose  the  lower  sign.  Hence 

(s  - Y J 


B-.l  £n 


(10.14) 
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which  can  be  rearranged  as  follows.  We  have  (s  + v ) (s  - v ) ~ » 

s s o 


hence 


3 = -1  iln 


s + Y 


= +1  £n 


s + Y 


(10.15) 


The  definitions  (10.14)  and  (10.15)  of  the  function  cos”^ (s/k^) 

are  not  those  given  in  most  books  on  mathematical  functions.  In  those 

books  the  square  root  function  is  usually  understood  to  have  a branch 

cut  from  -k  to  +k  , and  is  quite  different  from  the  function  y which 
o o s 

occurs  in  wave  applications.  In  particular,  if  the  cut  goes  from  -k 

J.  ° 

to  +k  , (s^  - k^)^  behaves  like  s both  as  s -*•  + “ and  s whereas 


Y behaves  like  s as  s but  as  - s when  s 

s 


With  the  definition 


cos  ^ (-j^.  = + i Jin 
o 


("?) 


(10.16) 


c:onsider  the  functions  ^(s) , P_(s)  defined  in  (10.11).  In  the  first 

place,  no  new  branch  cuts  are  introduced  by  the  logarithm.  A branch 

cut  would  be  needed  only  If  s + Y =0  were  possible  for  some  value  of  s, 

s 

anf  no  such  value  of  s exists.  Thus  the  only  singularities  are  the  branch 
points  at  s = ± k^.  Consider  the  function  P^(s)  near  the  point  s = + k^. 


writing  s = k + u where  u is  small.  We  have 
o 


P (k  + u) 
+ o 


iu^(2k  + u)^ 
o 


i_ 

2 /OA,  a.  2 


uM2k  + u) 

i 

o o 


1 11.  1 

lu^(2k  )^  u2(2k  )2 

o o 


+ 0(u2) 


(10.17) 
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as  u -»■  0,  so  that  P.  (s)  is  in  fact  single-valued  near  s « k . The  branch 
+ o 

point  singularity  of  the  cos“^  cancels  that  of  Y at  s * k , and  hence 

s o 

P^(s)  Is  analytic  In  R^.  Similarly  P (s)  Is  analytic  In  R , the  branch 
point  singularity  of  Y at  s = - k being  cancelled  by  that  of 

8 O 

cos~^(-  s/k^)  there. 


To  verify  (10.12)  we  have,  for  s In  D, 


P^(s)  + P_(s) 


- s + Y, 


£n(-  1)  = Y, 


since  with  (10.13)  )ln(-l)  = - lit. 

To  find  the  behavior  at  Infinity  we  note  that  as  |s|  -►  «o  In  R^ 


I,  K 1 

' 7 ••• 


If  the  approach  to  infinity  is  below  the  branch  cut  from  s = + k^, 


Ys~-«  ••• 


(10.18) 


(10.19) 


if  s goes  to  infinity  above  the  cut.  Since  P^(s)  has  ^ singularity  at 

s > + k It  does  not  matter  which  of  these  Is  used,  and  we  find 
o 


~ s In  0(s”^  iln  s) 


(10.20) 


These  properties  of  the  P^  functions  arise  In  a great  many  applications 


to  wave  problems.  Corresponding  results  for  the  convected  wave  square 
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root  ^ can  be  derived  by  making  a change  of  variables  to  transform  ^ 

S 8 

into  as  was  done  in  ^9. 

C.  Incon  tressible  Flow  Problems 

Incompressible  flow  results  follow  from  taking  the  limit  k^-^  0, 

which  has  the  unfortunate  effect  of  reducing  the  strip  of  analyticity  D 

to  a line  on  which  the  functions  are  continuous,  but  not  necessarily 

analytic.  The  branch  cuts  from  ± k^  also  join  up  to  form  a complete 

barrier  along,  say,  the  imaginary  axis.  To  avoid  possible  difficulties 

stemming  from  this  it  is  usual  to  work  with  finite  k^  and  then  let  k^  -►  0. 

This,  however,  makes  for  unnecessary  complications  in  much  of  the  work, 

and  it  is  useful  to  be  able  to  tackle  the  incompressible  problem  directly. 

To  this  end  we  imagine  .the  branch  cuts  as  starting  from  0 + iO  and  going 

to  infinity  in  Im  s > 0 and  from  O-iO  to  infinity  in  Im  s < 0. 

Then  the  limiting  form  of  the  function  is  real  and  positive 

on  the  whole  real  s-axis,  i.e.,  it  is  there  the  function  (s|.  We  shall 
1 

write  (s^)^  for  this  function  in  the  complex  plane  with  the  cuts  as 
indicated;  it  is  the  continuation  to  complex  s of  the  function  |s|  on  the 
real  axis,  and  can  also  be  defined  as  s(sgn  Res)  where  sgn  x * ± 1 for 
x"^  0.  Thus 

1 

Yg  -*■  (s^)^  * s (sgn  Res)  as  k^  -*■  0 (10.21) 

The  multiplicative  split  is 
1 1 1 

(s^)'^-sTs‘^  (10.22) 

1 “ 1 i 

where  s^  means  the  branch  of  which  behaves  like  s^  as  s -►  +<«>with  a 

cut  from  0 - Oi  in  the  lower  half-plane,  while  behaves  like  s^  as 

s » but  has  the  cut  from  0+01  in  the  upper  half  plane. 
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The  additive  split  of  (s^)^  can  be  found  by  taking  the  limit  of 

(10.11)  as  k -*-0.  Define  In^s  to  be  the  branch  of  In  s with  a cut  from 
o + 

0-01  In  the  lower  half-plane,  and  with  in^s  real  and  positive  when  s Is 
real  and  greater  than  1.  Define  £n_s  similarly  except  that  It  has  the 
cut  from  0 + 01  In  the  upper  half-plane.  Then  we  can  see  that 

In^s  - in^s  “0  If  Res  > 0 

In^s  - i!,n_s  ■ 2iTT  If  Res  < 0 , 
and  both  of  these  are  covered  by 


In^s  - f.n_8  * 2iTT  H(-Res) 

where  H(x)  Is  the  Heaviside  function  equal  to  1 or  0 according  as 
X > 0 or  X < 0.  Since  sgn  x * 2H(x)  - 1,  we  can  write  this  now  as 

In^s  - In  8 ■ Itt  - liT  sgn  Res 
and  multiplying  by  s and  using  (10.21)  gives 

1 

s In^s  - 8 In  8 ■ iirs  - iit(8^)^ 

Now  define 


P.(s)  - f - ^ 


and  then  It  follows  from  (10.23)  that 
1 


(8*)*  - P^(8)  + P_(8) 


(10.23) 


(10.24) 


(10.25) 


and  P_^(8)  Is  analytic  In  R_j_(lm  s > 0),  P_(8)  Is  analytic  in  R_(Im  s < 0) , 
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as  required.  Further,  by  careful  consideration  of  the  branches  It  can 
be  shovm  that  these  P^(s)  still  have  the  property 

P_^(-s)  - P_(s) 

enjoyed  by  the  function  (10.11)  for  j*  0. 

Use  of  these  functions  enables  Incompressible  problems  to  be 
solved  much  more  elegantly  than  by  taking  the  limit  as  k of  the 


(10. 


more  complicated  compressible  problems,  but  we  should  stress  that  because 
there  Is  no  strip  of  overlap  all  procedures  should  only  be  regarded  as 
formal,  and  the  results  should  be  verified  by  Independent  checks. 


and  Strip  of  Analyticity  D 
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